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ABSTRACT 

Until recently, infrared observations of exoplanetary atmospheres have typically been interpreted us- 
ing models that assumed solar elemental abundances. With the chemical composition fixed, attempts 
have been made to classify hot Jupiter atmospheres on the basis of stellar irradiation. However, 
recent observations have revealed deviations from predictions based on such classification schemes, 
and chemical compositions retrieved from some datascts have also indicated non-solar abundances. 
The data require a two-dimensional characterization scheme with dependence on both irradiation and 
chemistry. In hot hydrogen-dominated atmospheres, the C/0 ratio critically influences the relative 
concentrations of several spectroscopically dominant molecules. Between a C/0 of 0.5 (solar value) 
and 2, the H2O and CH4 abundances can vary by several orders of magnitude in the observable at- 
mosphere, and new hydrocarbon species such as HCN and C2H2 become prominent for C/0 > 1, 
while the CO abundance remains almost unchanged. Furthermore, a C/0 > 1 can preclude a strong 
thermal inversion due to TiO and VO in a hot Jupiter atmosphere, since TiO and VO are naturally 
under-abundant for C/0 > 1. We, therefore, suggest a new two-dimensional classification scheme 
for hydrogen-dominated exoplanetary atmospheres with irradiation (or temperature) and C/0 ratio 
as the two dimensions. We define four classes in this 2-D space (01, 02, CI, and C2) with distinct 
chemical, thermal and spectral properties. Based on the most recent observations, we characterize the 
thermal structures and C/0 ratios of six hot Jupitcrs (XO-lb, CoRoT-2b, WASP-14b, WASP-19b, 
WASP-33b, and WASP-12b) in the framework of our proposed 2D classification scheme. While the 
data for several systems in our sample are consistent with carbon-rich atmospheres (i.e. C/0 >1), 
new observations are required to conclusively constrain the C/0 ratios in the dayside and terminator 
regions of their atmospheres. We discuss how observations using existing and forthcoming facilities 
can constrain C/0 ratios in exoplanetary atmospheres. 

Subject headings: planetary systems — planets and satellites: general — planets and satellites: indi- 
vidual 



1. INTRODUCTION 

Recent advancements in observational methods are fa- 
cilitating detailed characterization of exoplanetary atmo- 
spheres. Observations of thermal emission have been re- 
ported for over two dozen transiting exoplanets (Seager 
& Doming 2010), and have been used to constrain at- 
mospheric temperature profiles (Burrows et al. 2007, 
Knutson et al. 2008, Madhusudhan & Seager 2009,2010, 
Christiansen et al. 2010) and chemical compositions 
(Fortney et al. 2005, Seager et al. 2005, Barman et 
al. 2007, Grillmair et al. 2008, Madhusudhan & Sea- 
ger 2009,2011, Swain et al. 2009, Stevenson et al. 2010, 
Madhusudhan et al. 2011) of a substantial number of 
giant exoplanets. Similar successes, albeit for a smaller 
set of planets, have also been reported in characteriz- 
ing exoplanetary atmospheres using transmission spec- 
troscopy (e.g. Charbonneau et al. 2002, Vidal Madjar 
et al. 2004, Tinetti et al. 2007, Swain et al. 2008, 2012; 
cf Gibson et al. 2011, Bueaulieu et al. 2008, Desert et 
al. 2009, Rcdfield et al. 2008, Madhusudhan & Seager 
2009, Sing et al. 2011, Bean et al. 2010, Crossfield et al. 
2012, Berta et al. 2012), thermal phase curves (Knutson 
et al. 2007, 2012; Showman et al. 2009; Cowan et al. 
2012), and direct imaging (e.g. Marois et al. 2008,2010; 
Bowler et al. 2010; Currie et al. 2011; Janson et al. 
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2011). The majority of exoplanets whose atmospheres 
have been observed, both with ground-based as well as 
space borne facilities, are highly irradiated giant planets 
('hot Jupiters') whose large sizes and high temperatures 
make their spectra particularly favorable to detect using 
existing infrared instruments. 

Major efforts have also been directed towards develop- 
ing a consistent theoretical framework within which to 
interpret observed spectra of exoplanetary atmospheres. 
Atmospheric spectra are a function of the temperature 
structure and the chemical composition of the atmo- 
sphere. Early studies typically assumed solar abundances 
and chemical equilibrium in models of exoplanetary at- 
mospheres and focused primarily on characterizing hot 
Jupiters based on their temperature structures (Seager 
& Sasselov 1998; Sudarsky et al. 2003; Hubeny et al. 
2003; Barman et al. 2005; Fortney et al. 2008). Hubeny 
et al. (2003) demonstrated that highly irradiated at- 
mospheres which can be hot enough to contain gaseous 
TiO and VO, can host thermal inversion^ in them due 
to strong absorption of incident visible light by TiO and 
VO. The presence or absence of a thermal inversion, for a 
fixed chemical composition, causes distinct spectral sig- 
natures which can be identified in observations of thermal 
emission from hot Jupiters. On this basis, Fortney et al. 

^ A thermal inversions is a region in the atmosphere where 
temperature increases outward. Sec e.g. Madhusudhan & Seager 
(2010). 
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(2008) proposed two classes of irradiated giant planets 
based on their incident irradiation. Planets in the hotter 
class were predicted to host thermal inversions in their 
atmospheres, caused by the presence of gaseous TiO and 
VO; whereas those in the cooler class were predicted to 
lack thermal inversions in their atmospheres due to a lack 
of TiO and VO; the boundary between the two classes 
was chosen somewhat arbitrarily at an incident flux of 
~ 10^ ergs s^^cm^^. 

Subsequent observations and theoretical efforts sug- 
gested modifications to the TiO/VO hypothesis. While 
some highly irradiated hot Jupiters do seem to host ther- 
mal inversions (Burrows et al. 2007; Knutson et al. 2008; 
Christiansen et al. 2010; Madhusudhan & Seager 2010), 
several counterexamples to the TiO /VO hypothesis have 
also been reported. For example, some of the least irradi- 
ated hot Jupiters are reported as hosting thermal inver- 
sions (Machalek et al. 2008,2009), whereas several highly 
irradiated hot Jupiters show no conclusive evidence for 
thermal inversions (Fressin et al. 2010; Madhusudhan et 
al. 2011a; Anderson et al. 2012; Blecic et al. 2012). 
Spiegel et al. (2009) suggested that TiO and VO, be- 
ing heavy molecules, may be depleted from the upper 
atmospheres due to gravitational settling coupled with 
condensation. Consequently, they suggested, TiO may 
be responsible for causing thermal inversions in only the 
most irradiated of hot Jupiters, substantially hotter than 
the pM-pL boundary of ~ 10^ ergs s~^cm~^ originally 
proposed by Fortney ct al. (2008). While the argument 
of Spiegel et al could explain the lack of thermal inver- 
sions in some hot Jupiters, it does not explain planets 
with low irradiation which are reported to host thermal 
inversions (e.g. Machalek et al. 2008). 

More recently, Knutson et al. (2010) reported an em- 
pirical correlation suggesting that the chromospheric ac- 
tivity of host stars may inhibit the formation of ther- 
mal inversions in hot Jupiter atmospheres by potentially 
destroying inversion causing absorbers in the planetary 
atmosphere. While several planets have been know to 
conform to the correlation, challenges remain. Most of 
the inferences of thermal inversion considered in the cor- 
relation arc based on models with fixed solar chemical 
abundances (e.g. Todorov et al. 2012), while it has been 
shown that the temperature structure is highly degen- 
erate with chemical composition (Madhusudhan & Sea- 
ger 2010). For example, observations of two hot Jupiters 
TrES-2 and TrES-4 which were originally thought to host 
thermal inversions in the Knutson et al. study were sub- 
sequently found to be consistent with models without 
thermal inversions (Madhusudhan & Seager 2010). A 
similar example of XO-lb is also demonstrated in the 
present study. Nevertheless, any classification scheme 
for irradiated giant planets would have to consider the 
role of stellar activity, as suggested by Knutson et al, 
in destroying any inversion causing species that may be 
possible for a given chemistry. 

In view of the above challenges in characterizing hot 
Jupiter atmospheres, recent observations indicate the re- 
quirement of a fundamentally new dimension for charac- 
terizing exoplanetary atmospheres ~ a dimension which 
can represent chemistry instead of assuming solar abun- 
dances in chemical equilibrium. A key implication from 
all the previous studies (Hubeny et al. 2003, Fortney 
et al. 2008, Spiegel et al. 2009, Knutson et al. 2010) 



is that if the irradiation is high enough and the star is 
not active, TiO and VO would be prevalent in the atmo- 
sphere. Such an implication arises from the fundamental 
assumption of solar abundances in the atmospheric mod- 
els on which the above studies based their conclusions. 
Solar abundances imply an inherently oxygen-rich envi- 
ronment, with a carbon-to-oxygen (C/0) ratio of 0.5, 
due to which TiO and VO are abundant at high temper- 
atures. However, recent observations suggest that exo- 
planetary atmospheres may also be carbon-rich, i.e. C/0 
> 1 (Madhusudhan et al. 2011a; Cowan et al. 2012; 
Swain et al. 2012). For a C/0 > 1, it has been shown 
that TiO and VO are naturally under-abundant even at 
very high temperatures (Helling & Lucas 2009; Mad- 
husudhan et al. 2011b), thus making any of the other 
depletion processes such as due to gravitational settling 
or stellar activity irrelevant. Furthermore, a high C/0 
ratio fundamentally alters the chemical composition of 
the atmosphere at high temperatures by depleting strong 
opacity sources such as H2O and introducing new sources 
such as CH4, HON, and C2H2 (Seager et al. 2005; Mad- 
husudhan et al. 2011b; Kopparapu et al. 2012; Moses 
et al. 2012). Consequently, the C/0 ratio can have a 
substantial infiuence on the spectroscopic signatures of 
an exoplanetary atmosphere by influencing both its tem- 
perature structure as well as its chemistry. 

Motivated by recent observations, in the present work 
we suggest a new two-dimensional classification scheme 
in which the incident irradiation and C/0 ratio form the 
two dimensions. In this 2-D phase space, we identify four 
classes of H2-dominated atmospheres, namely 01, 02, 
CI, C2, each with distinct chemical, thermal, and spec- 
troscopic properties. While two of the classes (01 and 
02) pertain to oxygen-rich (C/0 < 1) chemistry, analo- 
gous to the two classes proposed by Fortney et al. (2008), 
the remaining two classes (CI and C2) host carbon-rich 
(C/0 > 1) chemistry. We discuss in detail the molec- 
ular constitution of atmospheres in the different classes 
and their corresponding spectral signatures, and apply 
the classification scheme to six hot Jupiters (XO-lb, 
CoRoT-2b, WASP-14b, WASP-19b, WASP-12b, WASP- 
33b) with different irradiation levels. We identify a sam- 
ple of 20 transiting exoplanets which form ideal candi- 
dates for characterizing their atmospheric C/0 ratios, 
and populating the 2-D phase space, and for which ob- 
servations are already available in at least three chan- 
nels of Spitzer photometry. By combining the existing 
observations with new near-infrared observations from 
ground and space, detailed constraints on the chemistry 
of atmospheres, interiors, and formation environments, 
can be obtained for a sizable sample of cxtrasolar plan- 
ets. Wc find in the present work that non-solar C/0 ra- 
tios can potentially explain observations of hot Jupiters 
which have hitherto been reported as anomalous using 
solar abundance models. 

Recently it has become possible to estimate elemental 
abundance ratios, such as C/H, 0/H, and C/0 (Mad- 
husudhan et al. 2011a) thus facilitating the thermal and 
chemical characterization of hot Jupiters in the proposed 
2D classification system in the future. Such elemental 
abundances in planetary atmospheres also provide criti- 
cal constraints on the chemistry of their interiors and on 
their formation conditions. In the solar system, the ob- 
served abundances of volatile elements such as C, S, N, P, 
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and noble gases in Jupiter's atmosphere support its for- 
mation by core accretion of planetesimals (Atreya et al. 
1999; Atreya & Wong, 2005) and imply a much cooler ori- 
gin (at temperature T < 30 K) for the planetesimals than 
Jupiter's present orbit at 5 AU where T ^ 160 K (Owen. 
1999). One of the most important quantities govern- 
ing the chemistry in protoplanetary environments is the 
C/0 ratio (Oilman 1969; Larimer 1974; Lodders & Feg- 
ley 1999; Lodders 2004; Kuchncr & Seager 2005; Bond et 
al. 2010; Madhusudhan et al. 2011b; Oberg et al. 2011). 
The chemical composition of the planetesimals and gas 
accreted by the forming planet can be oxygen-rich or 
carbon- rich depending on whether C/0 < 1 or C/0 > 1, 
respectively, in the planet-forming regions. However, the 
0/H and C/0 ratios of solar system giant planets are 
poorly known since their low atmospheric temperatures 
(T < 125 K) cause H2O to condense down to the deeper 
layers of the atmospheres (P '-^ 10 - 100 bar) which arc 
not observable in spectra obtained remotely. Neverthe- 
less, nominal constraints on the water abundances have 
been inferred indirectly from CO abundances (Fegley and 
Prinn 1988, Lodders and Fegley 1994, Bezard et al 2002, 
Visscher & Moses 2011). 

The 0/H and C/0 ratios are easier to measure for at- 
mospheres of highly irradiated extrasolar giant planets 
than they are for solar-system giant planets. The vast 
majority of exoplanets known have atmospheric tempera- 
tures of T - 500 3000 K (e.g. Charbonneau et al. 2005, 
2009; Deming et al. 2005; Knutson et al. 2008; Campo 
et et al. 2009; Smith ct al. 2011). At these tempera- 
tures, the major carbon and oxygen bearing molecules 
are in gaseous form in the observable atmosphere and 
hence can be inferred from spectra (Burrows & Sharp 
1999; Lodders & Fegley 2002; Madhusudhan & Seager 
2011). The spectroscopically dominant species are typi- 
cally H2O, CH4, CO, and CO2, allowing estimations of 
C/H, 0/H, and C/0 ratios. The detectability of these 
species depends on the spectral resolution and bandpass 
observed in, and on the molecular abundances, which 
depend on the atmospheric temperature-pressure profile 
of a given planet. While a statistically significant con- 
straint (at 3-(t) on the C/0 ratio of an exoplanetary at- 
mosphere has been possible only recently (Madhusudhan 
et al. 2011a), nominal constraints on atmospheric C, O, 
and/or C/O have been reported for several transiting ex- 
oplanets in the past decade (e.g. Fortney et al. 2005; Sea- 
ger et al. 2005; Madhusudhan & Seager 2009,2010,2011; 
Desert et al. 2009; Swain et al. 2009). However, cur- 
rently available facilities in space and on ground can be 
strategically used for high-confidence estimates of C/0 
ratios of a large sample of transiting exoplanets, as will 
be shown in the present work. The James Webb Space 
Telescope (JWST) will allow estimations of abundance 
ratios for a larger set of elements. 

In what follows, we describe our atmospheric model- 
ing methods in Section[5] We introduce the C/0 ratio in 
Section[3]as a viable new dimension for atmospheric char- 
acterization of exoplanets. Here, we study in detail the 
influence of the C/0 ratio on the chemical and spectro- 
scopic properties of cxoplanctary atmospheres and make 
observable predictions. We introduce a two-dimensional 
classification scheme for H2-dominated atmospheres in 
Section |4l In Section [5l we report model fits to obser- 
vations of six hot Jupitcrs, which show preliminary evi- 



dence for non-solar C/0 ratios. In Section [HI wc discuss 
observational approaches that might be most beneficial 
for characterizing cxoplanctary atmospheres using exist- 
ing facilities, and we identify currently known transiting 
exoplanets that are optimal for constraining C/0 ratios. 
Finally, we summarize our findings in section [71 

2. METHODS 

In this work, we study the dependence of the chemi- 
cal, thermal, and spectroscopic properties of hydrogen- 
dominated exoplanetary atmospheres as a function of two 
key variables - the C/0 ratio and the incident irradiation. 
Our primary focus in on irradiated transiting exoplanets, 
mainly hot Jupitcrs but our results are also applicable to 
hot Ncptunes, and hydrogen-dominated super-EarthfO in 
similar temperature regimes. While we compute models 
of cloud- free atmospheres, our results on gas phase chem- 
istry and corresponding spectroscopic signatures are also 
applicable to cloudy atmospheres. Planets at wide or- 
bital separations are expected to be cloudy, and their 
spectra contain absorption features of species resulting 
from gas phase chemistry over the cloud tops (e.g. Mad- 
husudhan et al. 2011c). In what follows, we describe 
our models for computing the atmospheric chemistry and 
spectra. 

We use the atmospheric modeling and retrieval tech- 
nique for exoplanetary atmospheres developed by Mad- 
husudhan & Seager (2009), henceforth MS09, and Mad- 
husudhan et al. 2011a. In order to compute spectra, 
the model treats 1-D line-by-line radiative transfer in 
a plane-parallel atmosphere, with constraints of local 
thermodynamic equilibrium (LTE), hydrostatic equilib- 
rium, and global energy balance, and includes the major 
molecular and continuum opacity sources. Since we are 
concerned with H2 dominated atmospheres, our primary 
sources of opacity in the infrared include, depending on 
the temperature and C/0 ratio, H2O, CH4, CO, CO2, 
C2H2, HCN, NH3, and H2-H2 CIA absorption. We also 
include TiO and VO in regions of the atmosphere where 
the temperatures exceed the condensation temperatures 
of TiO and VO (Spiegel et al. 2009). Our molecular line 
data are from Freedman ct al. 2008, Freedman (per- 
sonal communication, 2009), Rothman et al. (2005), 
Karkoschka & Tomasko (2010), and Karkoschka (per- 
sonal communication, 2011), and Harris et al. (2008). 
We obtain the H2-H2 collision-induced opacities from 
Borysow et al. (1997) and Borysow (2002). Given the 
parametric temperature and molecular abundance pro- 
files, the model computes a spectrum for the required 
geometry. For the majority of the work, we compute 
emergent spectra of thermal emission from the dayside 
atmospheres as observed at secondary eclipse, though in 
some instances we also compute transmission spectra as 
are observed at primary eclipse. 

We use our spectral model in two configurations. 
While computing sample spectra for representative ther- 
mal and chemical profiles we use our model in a for- 
ward modeling sense. On the other hand, we also use 
our models to derive best-fit solutions from observations 

It is presently unknown whether atmospheres of super-Earths 
arc hydrogen-rich or are dominated by heavy molecules, like ter- 
restrial atmospheres; observations have suggested both possibilities 
(Bean et al. 2010; Croll et al. 2011; Miller-Ricci et al. 2011). 
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Fig. 1. — Effect of C/O ratio on the volume mixing ratios of major species in H2-dominated atmospheres in chemical equilibrium. The 
solid curves show mixing ratios, with respect to H2, of H2O, CO, CH4, CO2, C2H2, C2H4, and HCN, at representative conditions (T ^ 2000 
K, P ~ 1 bar) in observable atmospheres of irradiated giant planets. A sharp transition in the abundances of all the molecules, except 
CO, occurs at C/O = 1. The dotted vertical line marks C/O = 0.54, corresponding to solar elemental abundances. Compared to solar 
abundance mixing ratios, H2O is depleted and CH4 is enhanced by over two orders of magnitude each for C/O > 1. 



of thermal emission from several hot Jupiters. In the 
later case, we use our model as a retrieval method (see 
Madhusudhan et al. 2011a). In this mode, the goal 
is to formally retrieve molecular abundances and tem- 
perature structures from given datasets, as implemented 
in Madhusudhan et al. (2011a) to estimate posterior 
probability distributions for the model parameters. This 
approach involves no assumptions about the chemical 
abundancef|3 or layer-by-laycr radiative equilibrium, and 
has been successfully used to constrain non-equilibrium 
chemistry and non-solar abundances (e.g. MS09, 2011; 
Stevenson et al. 2010; Madhusudhan et al. 2011a), and 
temperature inversions (MS09 & 2010; Christiansen et 
al. 2010) in several exoplanetary atmospheres. In order 
to sample the model parameter space, we use a Markov 
chain Monte Carlo (MCMC) algorithm with a Gibbs 
sampler (see Madhusudhan & Seager, 2010; Madhusud- 
han et al. 2011a). Given the limited amount of observa- 
tions available for each system, our goal in general is not 
to determine a unique best-fit model, but instead rule 
out regions of parameter space that arc excluded by the 
data. 

We place constraints on the temperature structure and 
chemical compositions of six hot Jupiters with different 
levels of incident irradiation: XO-lb, CoRoT-2b, WASP- 

* This approach is different from other forward models which 
typically assume solar abundances for exoplanetary atmospheres 
(e.g. Seager et al. 2005; Burrows et al. 2008; Fortney et al. 
2008). Once our solutions are retrieved from the data, they arc 
tested for physical plausibility using chemical equilibrium and non- 
equilibrium calculations and any deviations are reported. 



14b, WASP-19b, WASP-12b, and WASP-33b. In each 
case, we explore models with and without thermal inver- 
sions and with different C/O ratios. For each system, we 
identify sample models with different C/O ratios and/or 
thermal inversions that can explain the data in a self- 
consistent manner. And based on the model fits to the 
data we discuss new observations that might help fur- 
ther constrain the C/O ratio and thermal structure of 
each planet. 

Wc also investigate in detail the observable effects of 
the C/O ratio on the major molecular species. Simi- 
lar studies have been reported in the past by Kuchner 
& Seager (2005), Seager et al. (2005), and Madhusud- 
han et al. (2011b). For given elemental abundances, 
with a prescribed C/O ratio, and thermodynamic condi- 
tions (P and T) we compute the concentrations of the 
major species of carbon and oxygen in chemical equilib- 
rium. For this purpose, we use the equilibrium chemistry 
code originally developed by Seager et al. (2000), and 
subsequently used in several recent works (Seager et al. 
2005; Miller-Ricci et al. 2009; Madhusudhan & Seager 
2011). The code calculates gas phase molecular mixing 
ratios for 172 molecules, resulting from abundances of 23 
atomic species, by minimizing the net Gibbs free energy 
of the System. The multi-dimensional Newton-Raphson 
method described in White et al. (1958) is used for the 
minimization. The free energies of the molecules are cal- 
culated using polynomial fits based on Sharp & Huebner 
(1990). 

We study the influence of C/O ratio on all the spec- 
troscopically dominant species containing C and O in 
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Fig. 2. — Dependance of molecular mixing ratios on C/O ratio and temperature (T). Each panel is analogous to Fig. [T] for a particular 
T, shown in the legend, and P = 1 bar. At low T {< 1200 K), H2O is the dominant oxygen carrier and CH4 is the dominant carbon carrier, 
and CO is least abundant of the three; CO mixing ratio is below 10"^'' for T < 600 K. At high T (> 1200 K), the mixing ratios of H2O 
and CH4, and other hydrocarbons, are critically dependent on the C/O ratio, as discussed in Scction l3.1l 



hydrogen-dominated atmospheres, in chemical equihb- 
rium. Depending on the C/0 ratio those species include 
H2O, CH4, CO, CO2, C2H2, and HCN. Chemical equi- 
librium is a valid assumption for the highly irradiated 
atmospheres we consider in this study (e.g. Moses et 
al. 2011). However, for cooler planets non-equilibrium 
chemistry in the form of vertical eddy mixing and photo- 
chemistry can cause significant changes in the molecular 
mixing ratios, and hence, affect their spectral signatures 
(Zahnle et al. 2009; Line et al. 2010; Madhusudhan & 
Seager 2011; Moses et al. 2011; Vischer & Moses 2011). 
We consider the effect of C/0 ratio on non-equilibrium 
chemistry of hot Jupiter atmospheres in a companion 
study (Moses et al. 2012). 

3. C/O RATIO AS A DIMENSION FOR ATMOSPHERIC 
CHARACTERIZATION 

In this section, we study the influence of the C/0 ra- 
tio on the atmospheric properties of cloud-free hydrogen- 
dominated atmospheres. Spectra of such atmospheres 
are dominated by features of oxygen and carbon based 
molecules, O and C being cosmically the most abun- 
dant elements after H and He. For a given pressure- 
temperature (P-T) profile, the C/0 ratio is the most 
critical parameter that governs the concentrations of the 
dominant O and C bearing species, and, therefore, af- 
fects their spectral features and their influences on ther- 
mal profiles. Several studies have discussed the effect 
of the C/O ratio on these various aspects of hydrogen- 



dominated atmospheres (Kuchner & Seager 2005; Seager 
et al. 2005; Helhng & Lucas 2009; Madhusudhan et al. 
2011b; Kopparapu et al. 2012). Here, we explore a wider 
range of temperatures and C/0 ratios and re-evaluate 
the effects of the C/0 ratio from the standpoint of key 
observables. 

3.1. Effect of C/0 on C and O Chemistry 

At high temperatures, two distinct regimes of atmo- 
spheric chemistry can be identified based on the C/0 
ratio. Fig. [T] shows the mixing ratios of major O and 
C based molecules as a function of the C/0 ratio, for 
nominal conditions in the observable atmospheres of hot 
Jupiters (T = 2000 K, P = 1 bar). A sharp transition in 
the abundances of all the molecules, except CO, occurs 
at C/0 = 1. For example, in going from the oxygen- rich 
solar abundances with C/0 = 0.5 to a carbon-rich value 
of 1.5, the H2O mixing ratio drops by about three or- 
ders of magnitude and the CH4 mixing ratio increases 
by over three orders of magnitude. In addition, sev- 
eral hydrocarbons besides CH4, such as C2H2, C2H4, 
and HCN, also become markedly abundant for C/0 > 
1, whereas CO2 becomes negligible, as shown in Fig. [TJ 
These extreme changes represent the two distinctly ob- 
servable differences between oxygen- rich (C/O < 1) and 
carbon-rich (C/0 > 1) atmospheres at high tempera- 
tures, which make hot Jupiters ideal laboratories for con- 
straining C/0 ratios. On the other hand, at these T and 
P the CO abundance for C/0 between 0.5 and 1.5 re- 
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mains nearly unchanged in comparison. Consequently, 
the CO/H2O ratio, which is nearly unity for C/0 = 0.5, 
can attain values of up to 10 - 10^ for higher C/0 ratios. 

The influence of C/O on the molecular abundances is 
temperature-dependent, as shown in Fig. [2] The mix- 
ing ratios of the major species H2O, CO, and CH4 in 
a hydrogen-dominated atmosphere are governed by the 
following summary reaction: 

CH4 + H20^CO + 3H2. (1) 

At a nominal pressure (P) of ^1 bar, reaction ([1]) fa- 
vors the conversion of CH4 to CO at high T (> 1300 K), 
and the reverse reaction is favored at lower T . Conse- 
quently, CH4 is the dominant carbon-bearing molecule at 
low T and CO is the dominant carbon-bearing molecule 
at high T. In an atmosphere with a solar abundance com- 
position (C/0 = 0.5) which has surplus oxygen, most of 
the oxygen is contained in H2O at low T. At high T, 
about half the available oxygen pairs up with the car- 
bon atoms in CO, and the remaining half is present in 
H2O. Therefore, H2O is abundant at all T for a solar 
abundance atmosphere. On the other hand, for C/0 
> 1, H2O is still the dominant oxygen carrier at low T 
(< 1200 K) where CO is minimal, but at high T where 
CO is dominant, there is little excess oxygen available to 
form H2O, contrary to the solar abundance case. The 
excess carbon available for the C/0 > 1 case is present 
in the form of abundant CH4 and higher hydrocarbons 
such as C2H2, C2H4, and HCN. 
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Fig. 3. — Mixing ratios for carbon-rich chemistry relative to so- 
lar abundance chemistry as a function of temperature. The solid 
(dashed) curves show mixing ratios of H2O, CH4, and CO for C/O 
= 1 (C/O = 2) relative to those obtained for C/O = 0.5. The en- 
hancement of CH4 and depletion of H2O are apparent for T > 1200 
K, and can be up to three orders of magnitude for C/O between 1 
and 2, and for high T. The absolute mixing ratios are shown in|4] 

The relative changes in the H2O and CH4 abundances 
between the oxygen-rich and carbon-rich regimes as a 
function of temperature are shown in Fig. [3l The figure 
shows enhancements in CH4 and depletion in H2O for 
C/0 = 1 and C/0 = 2, compared to C/0 = 0.5 (solar 
value); the absolute mixing ratios for the different C/0 
ratios are shown in Fig. S] For T < 1200 K at P 1 
bar, CO is less abundant, and CH4 is the dominant car- 
bon carrier and H2O the dominant oxygen carrier, for all 
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Fig. 4. — Mixing ratios of H2O, CH4, and CO for the cases shown 
in Fig. O Mixing ratios are shown for C/O ratios of 0.5 (solar), 1.0 
and 2.0. H2O mixing ratios as low as 10^® and CH4 mixing ratios 
as high as IQ-^^ can be attained for C/O = 1 and 2, for T > 2000 
K in the lower atmospheres (P ~ 1 bar). 

C/0 ratios. In such atmospheres, estimating the C/0 
ratio requires high precision, to within a factor of 2, esti- 
mation of H2O and CII4 abundances; currently molecular 
abundances for exoplanetary atmospheres are measured 
to only within a factor of 10 at best (Madhusudhan & 
Seager 2009; Madhusudhan et al. 2011a). On the other 
hand, for T > 1200 K at P ^ 1 bar, as is the case for 
the majority of transiting exoplanets known, the CII4 
enhancement and H2O depletion increase with tempera- 
ture, and can be over two (three) orders of magnitude for 
C/0 = 1 (C/0 = 2) and T > 2000. In absolute terms, 
as shown in Fig. [H for C/0 > 1 the CH4 volume mix- 
ing ratios can be as high as ^ 10~^ even at T > 2500 
K, contrary to values < 10^^ predicted by a solar C/0 
ratio. At the same temperatures, for C/0 > 1, the II2O 
mixing ratios can be < 10~^, while solar abundances 
predict - 5 x lO""*. These CH4 and H2O mixing ratios 
for C/0 > 1 at high T are in agreement with the ob- 
servational constraints for WASP- 12b reported by Mad- 
husudhan et al. (2011a). Additionally, other hydrocar- 
bons such as C2H2, C2II4, and HCN, which had not been 
considered by Madhusudhan et al. (2011a), are also ex- 
pected to be present in the atmosphere of WASP-12b, as 
has also been suggested by previous studies (Madhusud- 
han et al. 2011b; Kopparapu et al. 2012). In Section [531 
we demonstrate that model spectra which include these 
species are consistent with spectroscopic observations of 
WASP-12b. 

3.2. Effect of C/0 on Thermal Inversions 

The C/0 ratio places important constraints on the 
likelihood of thermal inversions in exoplanetary atmo- 
spheres, as discussed in detail in Madhusudhan et al. 
(2011b). Here, we briefly summarize those results. Ther- 
mal inversions are caused by the presence of chemical 
species in the observable atmosphere which are strong 
absorbers of incident UV and visible light. Tradition- 
ally. TiO and VO have been proposed to cause thermal 
inversions in highly irradiated hot Jupiter atmospheres 
(Hubeny et al. 2003; Fortncy et al. 2008, but cf Spiegel 
et al. 2009). However, Madhusudhan et al. (2011b) show 
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that for C/0 > 1, the TiO and VO abundances in chem- 
ical equihbrium are lower by about two orders of magni- 
tude compared to their abundances for a solar C/0 ratio 
which was assumed in previous models of thermal inver- 
sions (Burrows et al. 2008; Fortney et al. 2008; Spiegel et 
al. 2009). Non-equilibrium effects such as gravitational 
settling can further deplete the abundances (Spiegel et al. 
2009). Additionally, stellar activity might also effect the 
photochemical stability of inversion-causing compounds 
(Knutson et al. 2010). Consequently, if C/0 > 1, even 
the most highly irradiated atmospheres cannot host ther- 
mal inversions due to TiO/VO, which might explain the 
lack of a thermal inversion in WASP-12b (Madhusudhan 
et al. 2011a, b). Conversely, very highly irradiated atmo- 
spheres which do not host thermal inversions may likely 
be carbon-rich. 

C/0 ratios derived from observations can help con- 
strain the chemical nature of inversion-causing absorbers. 
Besides TiO and VO, other oxygen bearing molecules are 
also unlikely to cause thermal inversions in C-rich atmo- 
spheres. For C/0 > 1, almost all the oxygen is occupied 
by CO, which is the dominant O bearing molecule at high 
T, leaving little O for other oxygen-bearing molecules. 
As such, if thermal inversions are detected in high C/0 
atmospheres, they must be formed either by hydrocar- 
bons or by other oxygen-free compounds, e.g. sulfur 
species (Zahnle et al. 2009). Hydrocarbon-based pho- 
tochemical products are responsible for thermal inver- 
sions in giant planets in the solar system, which are at 
much lower temperatures (T ~ 50 ~ 150) compared to 
hot Jupiters (Chamberlain 1978; Yung & Demore 1999). 
However, whether such absorbers can sustain in highly 
irradiated hot Jupiters (T ~ 1000 — 3000) remains to be 
investigated in future studies. Thus, a simultaneous con- 
straint on the C/0 ratio and the presence of a thermal 
inversion from a given data set is essential to characterize 
the chemistry of inversion causing absorbers in irradiated 
giant planet atmospheres. 

3.3. Characteristic P and T 

The measurability of the C/0 ratio in an cxoplanctary 
atmosphere depends on the characteristic temperature of 
the atmosphere. As shown in Figs. [2] and [U the varia- 
tions in the H2O and CH4 mixing ratios with the C/0 
ratio are significant only for T > 1200/^. Above 1200 K, 
the differences in the H2O and CH4 abundances between 
the 0-rich and C-rich regimes increase with T, reaching 
factors as high as 10^ (see Fig. [3]). Consequently, spectra 
of hotter atmospheres are generally more likely to carry 
stronger signatures of their C/0 ratios, and would hence 
be more effective in distinguishing between low (< 1) and 
high (> 1) C/0 models. 

The chemistry in the observable atmosphere (at P ~ 
0.01 - 1 bar; SI of Madhusudhan et al. 2011a) is in- 
fluenced by the temperature and pressure in the lower 
layers of the atmosphere. Vertical mixing in the atmo- 
sphere transports species from hotter lower regions of the 
atmosphere to upper regions, resulting in nearly uniform 
concentrations in the observable atmosphere (Cooper & 
Showman 2007; Zahnle et al. 2009; Line et al. 2011; 
Madhusudhan & Seager 2011; Moses et al. 2011). There- 
fore, the pressure level from where the species can be 
transported up (called the 'quench' pressure, Pg), and the 
temperature {Tq) there, are the equilibrium conditions 



governing the observed concentrations of the species. 

While Pq and Tq can be computed for planet-specific 
models based on non-equilibrium chemistry (e.g. Mad- 
husudhan & Seager 2011; Moses et al. 2011), we adopt 
generic values for our current purpose. For irradiated gi- 
ant planet atmospheres, Pq typically lies in the 1-10 bar 
range (Cooper & Showman, 2007; Madhusudhan & Sea- 
ger 2011; Visscher & Moses 2011). We assume a nominal 
value of Pq = 1 bar; larger Pq requires stronger vertical 
mixing. As for Tq, the dayside P-T profiles of strongly 
irradiated atmospheres asymptote to isotherms at high 
optical depth, for P ^ 1 — 10 bar (Burrows et al. 2008; 
Hansen 2009; Guillot 2010). The corresponding tem- 
perature (Tq) is a function of the visible and infrared 
opacities, the Bond albedo, and day-night energy redis- 
tribution (Hansen 2008; Guillot et al. 2010; Heng et al. 
2011), and can be quite high (up to a factor of ^ 1.4) 
compared to the equilibrium temperature which is given 

by 

Teq = r.[(l - Ab){1 - U){Rl/2a^)]"\ (2) 

where, T-,, and P^ are the effective temperature and ra- 
dius of the host star, and a is the orbital separation. Ab 
is the Bond albedo and is the fraction of incident flux 
redistributed to the night side (Madhusudhan & Seager 
2009). Consequently, based on comparisons with pub- 
lished profiles of several irradiated atmosphere models 
(Burrows et al. 2008, Fortney et al. 2008, Madhusudhan 
& Seager 2009), we find that a lower limit on Tq can be 
approximated by: 

T^>T,[Rl/2a'Y'^ = T,^\f=o, (3) 

where / = stands ioY Ab = fr = 0. 

We note that our estimation of Tq is only nominal, 
since several factors of non-equilibrium chemistry, such 
as vertical eddy mixing and photochemistry, can play a 
role in determine the quench levels of the species and Tq 
can be different for different species (Moses et al. 2011). 

The above estimates for Pq and Tq can be used to 
identify systems that are most suitable for C/0 mea- 
surements. Irradiated atmospheres with Tq > 1200 K at 
Pq = 1 bar can be expected to contain measurable sig- 
natures of their C/0 ratios, as shown in Figs. [2] and |3l 
The signatures are enhanced with temperature, and are 
maximimum for T > 2000 K. 

3.4. Effect of C/0 on Atmospheric Spectra 

The chemistry in C-rich (C/0 > 1) and 0-rich (C/0 
< 1) atmospheres present distinct spectroscopic signa- 
tures that are observable with current and forthcoming 
instruments. Atmospheres of transiting planets can be 
observed in transit spectra as well as in thermal spectra 
at occultatioi0. A transit spectrum constrains the chem- 
istry at the limb of the planetary atmosphere, whereas 
a thermal spectrum constrains the chemistry as well as 
the temperature structure in the dayside atmosphere of 
the planet (Madhusudhan & Seager 2009). 

In this section, we briefiy illustrate the effect on the 
C/0 ratio on the spectra of highly irradiated hydrogen- 
dominated planets. Models for cooler atmospheres have 

^ The spectral signatures of thermal spectra of transiting planets 
should also apply to spectra of directly imaged planets at similar 
effective temperatures, though the pressure-temperature profiles 
would be different. 
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Fig. 5. — Model transit spectra of hydrogen-rich atmospheres 
with different C/O ratios and temperature structures. The plane- 
tary properties (radius and surface gravity) are adopted from those 
of HD 189733b. Spectra are shown for C/O ratios of 0.5, 1.0, and 
2.0, for two different pressure-temperature {P-T) profiles. The up- 
per (lower) panel corresponds to a P-T profile with a temperature 
of 2000 (2600) K in the lower atmosphere (P = 1 bar). In models 
with C/O = 0.5 (red spectra), the peaks correspond to absorption 
features of H2O and CO; CH4 is negligible. For C/O = 1 (gray 
spectra), on the other hand, the H2O features are substantially 
diminished, and mild features of CH4 are visible, along with CO 
features. For C/O = 2, If20 features arc almost non-existent., and 
strong CH4 features arc apparent, in addition to CO features, as 
discussed in Section 13.41 



been reported by Kuchncr & Seager (2005) and Mad- 
husudhan ct al. (2011b). For ease of illustration, we 
consider only four molecules (H2O, CO, CH4, and CO2) 
in the models discussed here; detailed models includ- 
ing additional opacities due to HCN and C2H2 that arc 
possible in C-rich atmospheres will be demonstrated in 
Section [SJ Model transit and thermal spectra of tran- 
siting hot Jupiters with different C/O ratios and tem- 
perature profiles are shown in Figs. [5] & |6l Here, we 
use the planetary and stellar parameters of the canonical 
hot jupiter HD 189733b, and vary the temperature pro- 
file (adapted from Madhusudhan & Seager 2009) to span 
different high-temperature regimes, which is equivalent 
to changing the orbital separation. 

We generate model spectra for two different temper- 
ature profiles, shown in the insets in Fig. |6l with dif- 
ferent quench temperatures. The quench temperature 
(Tg), discussed in Section [331 is the temperature of the 
lower atmosphere (at P ~ 1 bar), which typically governs 
the chemistry in the observable atmosphere. Spectra are 

2600 K, and for C/O 



shown for = 2000 K and 

ratios of 0.5, 1.0, and 2.0, for each T„. For both T, 
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models with C/O = 0.5 show strong H2O features near 
1.1 /im, 1.4 /im, 1.8 /xm, and 3 ^m, which are clearly 
distinguishable from the lack of the same features for 
models with the C/O > 1. On the other hand, models 
with C/O > 1 have enhanced CH4 abundances which 
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Fig. 6. — Model thermal emission spectra of hydrogen-rich at- 
mospheres with different C/O ratios and temperature structures. 
Each panel shows spectra for a range of C/O ratios for the pressure- 
temperature (P-T) profile shown in the inset. The red, blue and 
green spectra correspond to C/O ratios of 0.5, 1.0, and 2.0, respec- 
tively. The influence of C/O ratios on the spectroscopic features 
are discussed in Section 13.41 



cause greater absorption in the CH4 bands near 1.2 ^m, 
1.3 /im, 2.3 /im, 3.3 /xm, and 7.8 /xm. However, the 
strength of the CH4 features depend on the temperature 
and C/O ratio. For Tq = 2000, the methane features 
are only marginally detectable for a model with C/O = 
1, but become the dominant spectroscopic features for 
C/O = 2. On the other hand, for 2600, only the 

strongest CH4 bands, at 3.3 /im and 7.8 /^m, are promi- 
nent for C/O > 1. The reason for the relatively weaker 
CH4 features for T ~ 2600, even for C/O > 1, is that 
even though the CH4 enhancement over solar models is 
as high as for T ^ 2000, the absolute CH4 abundances 
are lower by a factor of ~ 10, as shown in Fig. [21 Finally, 
all the models show strong CO features at 2.3 /im and 
4.8 /im, since CO is a prominent molecule in both C-rich 
and 0-rich regimes. 

The spectroscopic features discussed above can also be 
observed in thermal emission spectra observed at sec- 
ondary eclipse. Model spectra of hot Jupiters in thermal 
emission are shown in Fig. [51 In the absence of thermal 
inversions, spectra of high temperature 0-rich models 
show strong absorption in the H2O bands and weak ab- 
sorption in CH4 bands, and vice- versa for C-rich models. 



C/0 ratios in cxoplanctary atmospheres 



9 



It must be noted that for thermal spectra stronger ab- 
sorption means lower flux and hence a weaker signal. 
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Fig. 7. — A two-dimensional classification scheme for hydrogen- 
dominated atmospheres (see Section [4]|. The left and right verti- 
cal axes show the incident irradiation at the substeUar point and 
the equilibrium temperature with no redistribution and no al bedo , 
which is representative of the quench temperature (see Section [3.3p . 
Four classes are shown: Ol, 02, CI, and C2. Atmospheres in Ol 
and 02 are O-rich (C/O < 1), and those in CI and C2 are C- 
rich (C/O > 1). The distinction between Ol (CI) and 02 (C2) 
is based on irradiation, or temperature. The major molecules in 
each class are shown, for a pressure of 1 bar. Atmospheres in Ol 
and 02 classes are abundant in H2O. CH4 is negligible in 02, but 
can become significant in Ol-class atmospheres depending on the 
temperature; the CO - CH4 transiting temperature is ~ 1200 K 
(see Figs. [2] & EI . Being extremely irradiated, atmospheres in the 
02 class are more likely to host thermal inversions due to gaseous 
TiO/VO than the Ol-class atmospheres; classes Ol and 02 are 
analogous to the pL and pM classes of Fortney et al. (2008). The 
gray area in the 02 class indicates the uncertainty in incident ir- 
radiation at which gaseous TiO/VO can remain aloft in the upper 
atmospheres (cf. Spiegel et al. 2009). Atmospheres in CI and C2 
classes are under- abundant in H2O compared to their O-rich coun- 
terparts at the same irradiation. In C2-class atmospheres, H2O 
is negligible (< 10"'^), where as CO, CH4, C2H2, and HCN are 
abundant. In Cl-class atmospheres, H2O is significant, but still 
less than their O-rich counterparts. In both CI and C2 classes, 
TiO and VO are naturally under-abundant even for very high ir- 
radiation levels and, hence, cannot cause thermal inversions (Mad- 
husudhan et al. 2011b). 



4. A TWO-DIMENSIONAL CLASSIFICATION SCHEME 

The distinct changes in molecular composition with 
the C/0 ratio and their influence on the observable 
spectroscopic properties make the C/0 ratio a viable 
dimension for characterization of hot exoplanetary at- 
mospheres. Consequently, we suggest a two-dimensional 
classification scheme for highly irradiated giant planet at- 
mospheres, in which incident irradiation and chemistry, 
represented by the C/0 ratio, are the two dimensions. 
Our approach is a generalization of the one-dimensional 
scheme proposed by Fortney et al. (2008) in which irra- 
diation is the only dimension and the composition was 



fixed to be solar (C/0 = 0.54). Our addition of the C/0 
ratio as an extra dimension is also motivated by the pos- 
sibility that non-solar C/0 ratios might explain several 
of the observations which would otherwise be deemed 
as anomalies, based on the 1-D hypothesis, as shown in 
Section [5] and Madhusudhan et al. (2011a). 

A schematic representation of the two-dimensional 
classification scheme is shown in Fig. [T] Hydrogen- 
dominated irradiated atmospheres can be classified into 
four classes in this 2-D phase space. Along the C/0 
axis, C/0 = 1 forms a natural boundary between C-rich 
and O-rich atmospheres with manifestly distinct chem- 
ical and spectroscopic characteristics, as shown in Sec- 
tion [3] Along the irradiation axis, atmospheres can be 
classified broadly in two classes based on the temper- 
atures where significant transitions occur in the chem- 
istry in each C/0 regime, similar to the 1-D classification 
scheme of Fortney et al. (2008) for O-rich atmospheres. 
Consequently, we identify the four classes as: 



• 01 

• 02 

• CI 

• C2 



C/0 < 1 and F < Fcnt 
C/0 < 1 and F > Fcnt 
C/0 > 1 and Tg < 1200 K 
C/0 > 1 and Tg > 1200 K 



In this notation, classes 01 and 02 correspond to 
oxygen-rich (C/0 < 1) atmospheres, and CI and C2 cor- 
respond to carbon-rich (C/0 > 1) atmospheres. In each 
C/0 regime, the suffix '1' represents the lower temper- 
ature (or irradiation) class and '2' represents the higher 
temperature class. As described in Section 13. 1[ for a 
given irradiation level C-rich atmospheres (in classes CI 
or C2) are depleted in II2O and enhanced in CH4, C2II2 
and HCN, compared to O-rich atmospheres at the same 
irradiation (in classes 01 or 02). Within CI and C2, 
the amounts of H2O depletion and CH4 enhancement, 
over solar values, increase with irradiation. In the O-rich 
classes, H2O and CO are abundant in both 01 and 02, 
except for the lesser irradiated atmospheres in 01 (with 
Tg < 1300 K), where CH4 can be a major carbon carrier 
and can almost entirely replace CO for Tg < 1000 K. 

The division between 01 and 02 along the irradia- 
tion axis concerns the possibility of thermal inversions 
in irradiated atmospheres. 01 and 02 are analogous to 
the pL and pM classes, respectively, of Fortney et al. 

(2008) . In the O-rich regime, very highly irradiated at- 
mospheres are likely to host gaseous TiO and VO which 
might cause thermal inversions. However, our division of 
01 and 02 at an irradiation of Fait = 10^ ergs/s/cm^, 
following Fortney et al., is only nominal. Spiegel et al 

(2009) showed that TiO and VO are subject to strong 
gravitational settling and hence may not be present aloft 
in the atmosphere in significant quantities for all but 
the most extremely irradiated systems. Consequently, 
a more accurate boundary might be at Fcrit 5 x 10^ 
ergs/s/cm^. This uncertainty in i^crit which is depicted in 
the gray area, could potentially be constrained by future 
constraints of thermal inversions in 02 class systems. 

The division between CI and C2 is motivated by the 
abundance of H2O and CO. In C-rich atmospheres, TiO 
and VO are naturally under-abundant, as discussed in 
Section 13.21 and hence atmospheres in neither CI nor 
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C2 can host thermal inversions due to TiO and VO. 
Consequently, at the current time, there is no motiva- 
tion to divide the high and low irradiation regimes on 
the basis of thermal-invcrsion-causing absorbers. How- 
ever, two regimes can be defined in C-rich atmospheres 
based on the effect of temperature on the relative abun- 
dances of H2O, CO, CH4. As shown in Section [5TT1 the 
depletion of H2O and enhancement of CH4 in C-rich at- 
mospheres over their 0-rich counterparts are substantial 
only for T > 1200 K (i.e. in C2 class). Atmospheres 
in the C2 class would be enriched in hydrocarbons such 
as CH4, C2H2, and HCN, and be manifestly depleted in 
H2O. As such, C2 class planets are readily identifiable 
by their weak H2O absorption bands, and strong hydro- 
carbon bands. For Tg < 1200 K, C-rich atmospheres 
have H2O depletion only by factors of a few compared to 
0-rich atmospheres with the same irradiation. At such 
temperatures, CI class planets may not be readily dis- 
tinguishable based on H2O abundances alone, unless re- 
ally precise estimates of the H2O abundances become 
available. In CI class objects CH4 is the most abundant 
molecule, followed by H2O and CO. Additionally, given 
their low temperatures CI class objects are most suscep- 
tible to non-equilibrium chemistry (Moses et al. 2011), 
by products of which may aid in their spectroscopic char- 
acterization. 

Classification of a planet in the present 2-D scheme re- 
quires that its irradiation and atmospheric C/0 ratio be 
known. While the incident flux can be obtained from the 
stellar and orbital parameters, determining the C/0 ra- 
tio requires infrared observations of the planetary atmo- 
sphere. As shown in blue in Fig. [51 nominal constraints 
on the C/0 ratios, based on at least six channels of 
broadband photometry and/or spectrophotometry, have 
been reported for only four transiting planets to date: 
HD 189733b (Madhusudhan & Scager, 2009; Swain ct al. 
2009), HD 209458b (Madhusudhan & Seager 2009,2010), 
GJ 436b (Stevenson ct al. 2010; Madhusudhan & Seager 
2011), and WASP-12b (Madhusudhan et al. 2011). Even 
so, only WASP-12b has a statistically robust constraint 
of C/0 > 1, and, hence, can be classifled as a C2-class 
system. The C/0 ratio of GJ 436b is currently uncertain 
owing to anomalously low CH4 abundances (Stevenson 
et al. 2010; Madhusudhan & Seager 2011) which could 
not be explained by non-equilibrium chemical models to 
date (Line et al. 2011; Madhusudhan & Seager 2011). 
Also shown in Fig. |8l in red, are tentative constraints for 
several hot Jupitcrs considered in the present study (see 
Section [5]) . 

The 2-D classification scheme presents a testable hy- 
pothesis for existing and upcoming observations. Since 
WASP-12b belongs to the C2 class, we expect it to be 
depicted in TiO and VO. This prediction is consistent 
with the lack of a noticeable thermal inversion in WASP- 
12b reported by Madhusudhan et al. (2011). We revisit 
WASP-12b in Section [5|). Nevertheless, future observa- 
tions of WASP- 12b in transit in the near- infrared (be- 
tween 0.7 fiui and 1.1 fim) could place constraints on 
the TiO / VO abundances in its atmosphere (Desert et al. 
2008). Similarly, the lack of noticeable thermal inversions 
in HD 189733b ( Burrows et al. 2008, Charbonneau et 
al. 2008, Madhusudhan & Seager 2009) and GJ 436b 
(Stevenson et al. 2010; Madhusudhan & Seager 2011) 
are also consistent with predictions for 01 systems. For 
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Fig. 8. — Two-dimensional characterization of transiting exoplan- 
ets for which constraints on C/O ratios arc available. Sec Fig. [71 
for description of the classification scheme. The arrows show the 
constraints on the C/O ratios. The blue arrows show constraints 
reported in previous studies (see Section |4]l, and the constraints in 
red show tentative C/O ratios inferred for systems in the present 
study. A question mark against a planet name indicates that its 
classification is only suggestive and requires new observations for 
confirmation (see Section [SJl. 

the hot Jupiter HD 209458b, the C/0 ratio is not weU 
constrained, spanning both 0-rich and C-rich regimes, 
depending on the models used to interpret the data (Sea- 
ger et al. 2005, Burrows et al. 2008, Madhusudhan & 
Seager 2009,2010). Consequently, more data is required 
to characterize HD 209458b. In the present study, we 
now consider six hot Jupitcrs over a wide range of inci- 
dent irradiation and attempt to characterize their atmo- 
spheres in the framework of our 2D classification scheme 
proposed here. 

5. APPLICATION TO SPECTRAL OBSERVATIONS OF HOT 
JUPITERS 

Recent observations of the dayside atmospheres of sev- 
eral hot Jupiters have proved difficult to explain in the 
context of traditional atmospheric models and classifica- 
tion schemes which assume solar abundances (i.e. C/0 
= 0.5). On one hand, several highly irradiated hot 
Jupiter atmospheres which were previously predicted to 
host thermal inversions due to TiO and VO showed no 
conclusive evidence of thermal inversions in their spectra, 
e.g. WASP-12b (Madhusudhan et al. 2011), WASP-19b 
(Anderson et al. 2012), WASP-14b (Blecic et al. 2012). 
On the other hand, in at least one case a thermal in- 
version was inferred for one of the least irradiated hot 
Jupiters, XO-lb (Machalek et al. 2008). In yet another 
case, no models have thus far been able to explain the 
dayside emission observed for the hot Jupiter CoRoT-2b 
(Deming et al. 2011). 

In the present section, we demonstrate that all the 
above mentioned observations can be consistently ex- 
plained in the framework of the 2-D classification scheme 
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TABLE 1 

Chemical compositions of model spectra 
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^ The two rows for each system correspond to the two model spectra for each system. "ORP" and "CRP" stand 
for oxygen-rich planet and carbon-rich planet, respectively. An ORP is defined as a planet with an atmospheric 
C/O < 1, and a CRP is defined as a planet with an atmospheric C/O > 1 (Madhusudhan et al. 2011b). 



proposed in section |4l using models that do not restrict 
the atmospheric compositions to solar C/O ratios. We 
present model fits to observations of six hot Jupiters 
spanning a wide range in stellar irradiation: XO-lb, 
CoRoT-2b, X0-2b, WASP-14b, WASP-19b, WASP-12b, 
and WASP-33b. For each planet, we consider models 
with 0-rich (i.e. C/O < 1) as well as C-rich (i.e. C/O 
> 1) compositions, and discuss the degree of fit in each 
case. Depending on the degree of fit, we attempt to char- 
acterize each planet based on its incident irradiation and 
possible C/O ratio. Two sample spectra, one 0-rich and 
one C-rich, fitting the data are shown for each planet. 
The corresponding chemical compositions are shown in 
Table m 

5.1. XO-lh 

The hot Jupiter XO-lb (McCullough et al. 2006) has 
presented one of the extreme counterexamples to the 
TiO/VO hypothesis to date. The planet orbits a G dwarf 
on a 4 day orbit, and is one of the least irradiated hot 
Jupiters whose atmospheres have been observed to date, 
with an equilibrium temperature of ^ 1400 K. As such, 
Fortney et al. (2008) predicted that the daysidc atmo- 
sphere of XO-lb would not be able to host a thermal 
inversion since it would not be hot enough to maintain 
gaseous TiO or VO in the upper atmosphere. Subse- 
quently, Machalek et al. (2008) reported observations 
of thermal emission from XO-lb in the four channels of 
Spitzer ISlkC photometry (at 3.6/im, 4.5/.im, 5.8/.im, and 
8/xm). Interpreting the data using models of Burrows et 
al. (2008), Machalek et al. (2008) reported a thermal in- 
version in the dayside atmosphere of XO-lb. Given the 
low irradiation received at its substellar point, the infer- 
ence of a thermal inversion in XO-lb presented a stark 
violation of the TiO/VO hypothesis. 
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Fig. 9. — Observations and models of thermal emission from the 
dayside atmosphere of XO-lb. The blue circles with error bars 
in the main panel show Spitzer IRAC observations from Machalek 
et al. (2008), and the solid curves show two model spectra. The 
green curve shows a model with oxygen-rich (C/O = 0.5) compo- 
sition, and the red curves shows a carbon-rich model with C/O = 
1. The carbon-rich model provides a better fit to the data. The 
oxygen-rich model undcr-predicts the flux at 5.8 fim. The band- 
pass integrated model points are shown in the colored circles. The 
corresponding pressure-temperature (P-T) profiles, which do not 
have thermal inversions, are shown in the inset. The black curves 
at the bottom of the figure show seven photometric bandpasses: J, 
H, Ks, and Spitzer IRAC bands at 3.6 ^m, 4.5 fim, 5.8 fim, and 8 
fim. The dotted lines show three blackbody spectra of the planet 
at 1000 K, 1250 K, and 1500 K. 

One potential explanation to the apparent contradic- 
tion was proposed by Knutson ct al. (2010) who found 
a correlation between inferences of thermal inversions 
in hot Jupiters and the chomospheric activities of their 



12 



host stars. They concluded that planets orbiting chromo- 
spherically active stars did not show evidence for thermal 
inversions in their atmospheres, and those orbiting quiet 
stars showed evidence for thermal inversions. Given the 
low chromospheric activity of its host star, the thermal 
inversion in XOlb appeared to be consistent with the em- 
pirical hypothesis of Knutson et al. However, it is still 
unknown to date as to what molecules in the atmosphere 
could actually cause a thermal inversion in XO-lb that is 
strong enough to match the observations of Machalek et 
al. (2008), since gaseous TiO and VO would be clearly 
underabundant in such a cool atmosphere (Spiegel et al. 
2009). In order to create the thermal inversion in their 
models, Machalek et al. had to use an artificial opacity 
source in visible wavelengths with a parametric opacity, 
wavelength range, and location in the atmosphere (Bur- 
rows et al. 2008). Whether such a source is physically 
plausible in XO-lb is yet to be demonstrated. 

In the present work, we report a consistent explana- 
tion to the observations of Machalek et al. (2008). The 
observations and two model spectra are shown in Fig. |9] 
and the corresponding chemical compositions are shown 
in Table [T] The four IRAC observations are shown in 
blue. Firstly, we confirm the findings of Machalek et al. 
that a solar abundance model cannot fit the data with- 
out a thermal inversion, as shown in the green curve in 
Fig. [9l The model fits the data at 3.6 /im, 4.5 fim, and 
8 fiia to within the 1-a error bars, but under-predicts 
the 5.8 fim fiux by 3-(t. Our non-inverted solar abun- 
dance model provides a slightly better fit to the data 
than that of Machalek et al. (2008), while still grossly 
under-predicting the flux at 5.8 //m, in agreement with 
Machalek et al. Furthermore, the temperature struc- 
ture required for a best-fit non-inverted solar abundance 
model as shown in Fig.|9]is nearly an isotherm, implying 
that the spectrum is equivalent to a blackbody, as shown 
by the dotted lines, making the 0-rich composition ir- 
relevant. The requirement of a thermal inversion with 
a solar abundance model thus rests largely on the fiux 
observed in the 5.8 /iui IRAC bandpass which has strong 
spectral features of II2O. In a solar composition model, 
II2O is the dominant oxygen bearing species at the tem- 
peratures of XO-lb. If there is no thermal inversion in 
the atmosphere, in which case temperature monotoni- 
cally decreases outwards, a low fiux is observed in the 
5.8 fim IRAC band due to strong water absorption. The 
contrary observation of the high 5.8 ^m flux in the data 
implies, in the context of a solar abundance model, that 
it must be due H2O in emission rather than absorption. 
This is possible in the presence of a thermal inversiorQ. 

We report, however, that the observations of XO-lb 
can be explained without the requirement of an ad hoc 
thermal inversion if the atmosphere of XO-lb is carbon- 
rich (C/O > 1). A model spectra of XO-lb with a C-rich 
composition and a non-inverted temperature profile, is 
shown in red in Fig. [HI The model fits all the observa- 
tions to within the Icr uncertainties, on average. The key 
to explaining the high flux in the 5.8/.im channel, simul- 
taneously with the other photometric points, without a 
thermal inversion lies in the H2O and CII4 abundances 
expected in a C-rich atmosphere (see e.g. Madhusudhan 

® A detailed description of such inferences is presented in Mad- 
husudhan & Seager, 2010). 



et al. 2011). At the temperatures of XO-lb, C/O > 1 
causes 10 — 100 times depletion of II2O and ^ 10 — 100 
times enhancement of CH4 compared to a solar abun- 
dance model (with C/O — 0.5). Additionally, other hy- 
drocarbons such as HCN and C2II2 arc also expected to 
be abundant. The low II2O abundance causes low ab- 
sorption in the 5.8^m channel, leading to a high thermal 
flux emitted from the deeper layers of the atmosphere. 
On the other hand, the low fluxes in the 3.6 fj,m and 8 
^m channels are explained by enhanced CII4, HCN, and 
C2H2 absorption in both channels. The 4.5 ^m flux in 
the C-rich model is explained by absorption due to CO 
and HCN. Therefore, we classify XO-lb as a C2 class 
planet in the framework of the 2-D classification scheme 
presented in section |4l 

Future observations in the near-infrared with existing 
facilities will be required to confirm the C2 classifica- 
tion of XO-lb. Given that our current inference is based 
solely on the one photometric point at 5.8 /xm, new ob- 
servations are required to confirm the potentially C-rich 
nature of XO-lb. As shown in Fig. [9l the high C/O 
models fitting the IRAC data predict substantially higher 
fiuxes in the 2-3 /im range compared to a solar compo- 
sition model. This difference is largely due to different 
continuum emission between the two models; the high 
C/O solutions have higher temperatures in the low at- 
mosphere (-P ^ 1 bar) which are required to fit the high 
5.8 /im IRAC point. Such differences can be used to 
constrain the models using ground-based high precision 
photometry (e.g. Lopez-Morales, 2010; Croll et al. 2011) 
and spectrophotometry (Bean et al. 2010; Mandell et 
al. 2011; Crossfield et al. 2011), and space-borne spec- 
troscopy with HST in the near term, and with JWST 
in the long term. Note that precisions better that 10^* 
would be required to distinquish between the C-rich and 
0-rich models. 

Care must be taken, however, in the interpretation 
of the C/O of XO-lb using transmission spectroscopy. 
Given the low irradiation levels received by the planet, a 
transmission spectrum of XO-lb is likely to probe even 
cooler regions of the atmosphere. It is possible that 
the temperatures in such regions can be well below the 
1200 K minimum required for high- confidence estimates 
of C/O ratios (see Section 3.1). Consequently, spectra of 
the limb of the atmosphere may show absorption features 
of H2O even if the C/O > 1. In such cases, establishing 
the C/O ratio accurately would require estimating the 
H2O and CH4 abundances to precisions within factors 
of a few. Such precisions would require extremely pre- 
cise ( 10~^) spectroscopic observations which may not 
be feasible with current instruments. At present, one 
could attempt to derive such molecular abundances us- 
ing the HST NICMOS spectrum of XO-lb reported by 
Tinctti et al. (2010), whose reported constraints agree 
with a potential C-rich atmosphere. However, we refrain 
from interpreting their data in the present work since 
their data reduction is currently contested (Gibson et al. 
2010). 

5.2. CoRoT-2h 

The hot Jupiter CoRoT-2b (Alonso et al. 2008) 
presents yet another challenge for atmospheric charac- 
terization of hot Jupitcrs using traditional models and 
classification schemes. The planet orbits a late G dwarf 
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Fig. 10. — Observations and model spectra of thermal emission 
for CoRoT-2b. The four observations are shown in blue circles 
with error bars at wavelengths 2.15 fim (Alonso ct al. 2010), 3.6 
(Doming ct al. 2011), 4.5 fim, and 8 fj,m (Gillon ct al. 2010 
& Doming ct al. 2011). The green curve in the main panel shows 
a model with an oxygen-rich (C/O j 1) atmosphere, and the red 
curve shows a model with a carbon-rich composition (C/O > 1). 
The corresponding pressure-temperature profiles are shown in the 
inset. The temperature profiles do not have thermal inversions. 
The red and green colored circles show the corresponding models 
integrated in different photometric bandpasses, the filter functions 
for which are shown in black at the bottom. The dotted curves 
show three blackbody spectra of the planet at 1000 K, 1500 K, and 
2100 K, for reference. 

on a 1.7 day orbit and is one of the highly irradiated hot 
Jupiters predicted by to host a thermal inversion in its 
dayside atmosphere (Fortney et al. 2008). However, the 
host star is known to be very chromospherically active, 
which according to the Knutson et al (2010) hypothesis 
would suggest that there should not be a thermal inver- 
sion in its dayside atmosphere. Between the two compet- 
ing predictions, however, lies the uncomfortable fact that 
neither models with thermal inversions nor models with- 
out thermal inversions have been able to explain the ob- 
served spectrum of CoRoT-2b (Gillon et al. 2010, Dem- 
ing et al. 2011). Observations of dayside thermal emis- 
sion from CoRoT-2b showed an anomalously low flux in 
the S^m Spitzer IRAC channel, and high flux in the 4.5 
fim channel, both of which could not be explained by 
previous atmospheric models (Alonso et al. 2010, Gillon 
et al. 2010, Deming et al. 2011, Guillot et al. 2011). 
As an alternative explanation, Deming et al. (2011) re- 
quired either excess emission from tidal mass loss that is 
too high to be realistic or the presence of an unknown 
absorber that selectively absorbs at wavelengths below 5 
/xm. 

In the present work, we are able to explain all the ex- 
isting observations of CoRoT-2b. Figure [TU] shows the 
observations and two model spectra fitting the data, and 
the corresponding chemical compositions are shown in 
Table [T] The data include broadband photometric detec- 
tions of thermal emission from CoRoT-2b in four infrared 
channels: a ground-based Ks band at 2.1 /xm (Alonso et 
al. 2010), and three Spitzer IRAC bandpasses at 3.6 /xm 
(Deming et al. 2011), 4.5 /xm and 8 /xm (Gillon et al. 
2010; Deming et al. 2011). The models shown in Fig- 
ure [10] include an 0-rich model (C/O < 1) and a carbon- 



rich (C/O > 1). The corresponding temperature profiles 
are shown in the inset of Figure [TOl As shown in the 
P-T profiles, the observations rule out a thermal inver- 
sion in the dayside atmosphere of CoRoT-2b. The low 
flux in the 8 /xm IRAC channel, thought to be anomalous 
by previous studies, can be naturally explained by CH4, 
HCN, and C2II2 absorption in the carbon-rich model and 
by H2O absorption in the oxygen- rich model. The same 
opacity sources also explain the low flux in the 3.6 /xm 
channel. 

Existing observations of CoRoT-2b are less constrain- 
ing of its C/O ratio, implying that it could be classified as 
01 or C2 according to the 2D classification scheme pro- 
posed in section 21 although the C2 class is marginally 
favoured. We find that all the observations can be ex- 
plained to within the 1-cr uncertainties on average by a 
wide range of models with carbon-rich (C/O > 1) and 
oxygen-rich (C/O < 1) compositions. The Spitzer IRAC 
5.8 /xm flux that implied a substantial depletion of H2O in 
XO-lb is unavailable for CoRoT-2b, leading to a poorer 
constraint on the II2O abundance in the latter. Thus, it 
is expected that both II2 0-rich and Il20-poor solutions 
are possible for CoRoT-2b, implying low C/O (< 1) and 
high C/O (> 1) ratios, respectively. This degeneracy 
occurs because II2O, CII4, and HCN, all have absorp- 
tion features in the 3.6 and 8 /xm IRAC channels. A 5.8 
/im observation could have broken the degeneracy with 
an independent constraint on the H2O abundance; the 
IRAC 5.8 /im bandpass does not contain significant CH4 
features. However, since the C-rich model provides a 
marginally better fit to the data, especially at 8 /xm, we 
tentatively classify CoRoT-2b as a C2 class planet in the 
framework of our 2D classification scheme proposed in 
section IH 

New observations of CoRoT-2b would be required to 
conclusively constrain its C/O ratio. As shown in Fig- 
ure [ini the two spectral models differ substantially in 
their absorption features. In particular, observations us- 
ing the HST WFC3 instrument in the 1.1-1.7 /xm band- 
pass would be able to conclusively determine the H2O 
abundance based on the strong H2O feature in the O- 
rich model around 1.4 fim. In addition, ground-based 
observations with higher precision that that currently 
available in the Ks band would be able to better con- 
strain the lower atmosphere temperature profile of the 
atmosphere, which in turn could differentiate between 
the different models. Another alternative would be to use 
transmission spectroscopy and photometry from ground 
and space to constrain the atmospheric composition of 
the day-night terminator of the planet. 

5.3. WASP-Ub 

WASP-14b is a transiting hot Jupiter orbiting an F5 
main-sequence star at an orbital separation of 0.036 AU 
(Joshi et al. 2009). With an incident irradiation of 2.5 
X 10^ ergs/cm^/s, leading to a zero-albedo equilibrium 
temperature of ^ 2200K, WASP-14b belongs to the class 
of highly irradiated hot Jupiters that are expected to 
host strong thermal inversions based on the TiO /VO hy- 
pothesis of Fortney et al. (2008). However, using three 
channels of SpitzerlRAC photometry, at 3.6 /xm, 4.5 /xm, 
and 5.8 fim, Blecic et al. (2012) reported the lack of a 
thermal inversion in the dayside atmosphere of WASP- 
14b. The absence of a thermal inversion in WASP-14b 



14 



5 F 



C/O = 0.1 i 




X (/irn) 

Fig. 11. — Observations and model spectra of thermal emission 
from the dayside atmosphere of WASP-14b. The data arc shown 
in blue circles with error bars obtained in the Spitzer IRAC band- 
passes at 3.6 fim, 4.5 fim, and 8 fj,m (Blecic ct al. 2012). The solid 
green and red curves show two model spectra with different C/O 
ratios. Both models provide reasonable fits to the data, though 
the oxygen-rich model provides a marginally better fit. The tem- 
perature profiles corresponding to the two models are shown in the 
inset. A thermal inversion is ruled out by the data. The Spitzer 
bandpasses are shown in the solid black curves at the bottom of the 
plot. The dotted gray curves show the planet blackbody spectra 
at 1600 K, 2200 K, and 2600K. 

was evident primarily from a low brightness temperature 
of ~ 1600K in the 8 ^m channel compared to temper- 
atures above 2200 K observed in the lower wavelength 
channels at 3.6 /xm and 4.5 /im, indicating temperature 
decreasing outward in the atmosphere. 

While the lack of a thermal inversion in WASP-14b 
may contradict the prediction of Fortncy ct al. (2008), 
it is also possible that condensation and gravitational 
settling of TiO/VO (Spiegel et al. 2009) might preclude 
the formation of a thermal inversion. It is unclear if pho- 
tospheric activity of the host star may have caused the 
disappearance of TiO and VO. Knutson et al. (2010) 
find that host star WASP-14 is not very active imply- 
ing that the planet should have had a thermal inversion. 
But Knutson et al. (2010) also highlight the difficulty 
in accurately establishing the Ca II H and K indices for 
F stars, such as WASP-14b. Alternately, as shown in 
Madhusudhan ct al. (2011b) and in the present work, it 
may also be possible that the atmosphere of WASP-14b 
is carbon-rich (i.e. C/O > 1) in which case TiO and VO 
are naturally low in abundance. 

We confirm the conclusions of Blecic et al. (2012) that 
the limited data available cannot conclusively constrain 
the C/O ratio of the dayside atmosphere of WASP-14b. 
Two model spectra fitting the observations arc shown in 
Fig. (TTJ and Table [T] shows the model compositions. In 
the 0-rich model, the dominant absorption in the 8-^m 
and 3.6- fim IRAC channels is caused by II2O, whereas 
that in the C-rich model is caused by HCN, CH4, and 
C2H2. Both models contain absorption features due to 
CO in the 4.5 fim bandpass. The temperature profiles for 
the two models are shown in the inset in Fig. 1111 Since 
the data rule out a thermal inversion in the atmosphere 
regardless of the C/O ratio, both models have similar 
lower atmospheric temperatures. 

Given that both C-rich and 0-rich models explain the 
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Fig. 12. — Observations and model spectra of the dayside at- 
mosphere of WASP-19b. The observations, shown in blue circles 
with error bars, comprise of three ground-based photometric data 
at 0.9 fim (Burton et al. 2012), 1.6 fim (Anderson et al. 2010) and 
2.1 fim (Gibson et al. 2010), and four photometric observations in 
the Spitzer IRAC bandpasses at 3.6 ^m, 4.5 fim, 5.8 ^m, and 8 
fim (Anderson et al. 2012). The solid red and green curves in the 
main panel show two model spectra with different C/O ratios as 
described in the legend. The data marginally favor a carbon-rich 
model (C/O > 1), though an oxygen-rich model (C/O < 1) cannot 
be conclusively ruled out. The corresponding temperature profiles 
are shown in the inset. A strong thermal inversion is ruled out 
by the data. The dotted lines show three blackbody spectra at 
1800 K, 2500 K, and 2900 K. The black curves at the bottom of 
the figure show the photometric bandpasses corresponding to the 
data. 

data almost equally well, WASP- 14b could be classified 
either as an 01 or a C2 system. Note that despite the 
high irradiation, WASP-14b cannot be classified as an 02 
system because of the lack of a strong thermal inversion 
and hence falls in the gray area between the 01 and 02 
classes shown in Fig. [71 However, since an oxygen-rich 
model (in green) provides a marginally better fit to the 
8- fim IRAC point than the carbon- rich model (in red), 
one may tentatively classify WASP-14b as an uppcr-01 
class planet, as shown in Fig. |8l 

New observations are needed to establish the C/O ratio 
of WASP-14b. The degeneracy between the two models 
arises from the lack of a spectral measurement in the 
IRAC 5.8 fim channel. In principle, a 5.8-fim IRAC data 
point, if available, could have easily removed the degen- 
eracy. Given that the 5.8- fim IRAC channel is no longer 
available on Spitzer, an alternative way to constrain the 
water abundance with existing facilities would be to ob- 
serve the II2O bands in the HST WFC3 bandpass be- 
tween 1.1 and 1.7 ^m. 

Future confirmation of the C/O ratio of WASP-14b 
would lead to one of two interesting conclusions. If the 
C/O is found to be less than unity, based on i/ST obser- 
vations, it will be a stringent constraint on the irradiation 
evidence for gravitational settling of TiO/VO according 
to the the Spiegel et al (2009), and will place a lower- 
limit on the irradiation level that defines the gray region 
between the 01 and 02 classes in Figs. [7] & [H Alter- 
nately, if the C/O is found to be > 1, WASP-14b wiU 
join WASP-12b in the new class of CRPs or carbon-rich 
planets (Madhusudhan et al. 2011b) . 

5.4. WASP-19b 
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The transiting hot Jupiter WASP-19b orbits a G8 
main-sequence star at a separation of 0.017 AU (Hebb ct 
al. 2010). At an irradiation level of 4 x 10^ ergs/cm^/s, 
it is one of the very highly irradiated hot Jupiters with a 
zero-albedo equilibrium temperature of 2400 K. Conse- 
quently, it is also expected to host a strong thermal in- 
version due to TiO / VO if the atmosphere is in chemical 
equilibrium and the composition is oxygen-rich as pre- 
dicted by Fortney et al. (2008). However, using six chan- 
nels of broadband photometry obtained using Spitzer and 
ground-based facilities, Anderson et al. (2012) report 
the non-detection of a thermal inversion in the dayside 
atmosphere of WASP-19b. Their inference was based 
on observations of thermal emission in four channels of 
SpitzerlRAC photometry at 3.6 ^m, 4.5 ^m, 5.8 /im, and 
8 /im, and two ground-based photometric measurements 
at 1.6 /im (Anderson et al. 2010) and 2.1 fim (Gibson et 
al. 2010). The inference of a non- inverted profile is evi- 
dent from the higher brightness temperatures observed in 
the 1.6-3.6 bands compared to the longer wavelength 
bands, suggesting a temperature decreasing outward in 
the atmosphere. 

Even though the lack of a noticeable thermal inver- 
sion in WASP-19b defies the predictions of the TiO/VO 
hypothesis of Fortney et al. (2008), alternate explana- 
tions exist. It might be unlikely that gravitational set- 
tling (Spiegel et al. 2009) alone can deplete all the TiO 
and VO in this highly irradiated atmosphere. However, 
Anderson et al. (2012) note that the planet orbits a chro- 
mosphcrically active star, in which case, according to the 
hypothesis of Knutson et al. (2010) the high UV flux in- 
cident on the planetary atmosphere could potentially dis- 
sociate the inversion causing molecules. An alternative 
explanation, as suggested by Madhusudhan et al. (2011), 
could be that the atmosphere is carbon-rich in which case 
TiO and VO would be naturally under-abundant and 
hence the lack of a thermal inversion. Therefore, WASP- 
19b could belong to either class 01/02 or C2, depending 
on the C/0 ratio. 

Current observations of WASP-19b are inconclusive 
about its C/0 ratio but nominally favor a C/0 > 1. The 
data and two model spectra of WASP- 19b are shown in 
Fig. [12] (sec Tabic [T] for the chemical compositions of the 
models). As demonstrated by Anderson et al. (2012), 
the majority of available data arc fit equally well by O- 
rich as well as C-rich models. The large error bar on 
the flux in the 5.8 /tm Spitzer IRAC channel makes it 
difficult to constrain the H2O abundance which is one of 
the major discriminator between the 0-rich and C-rich 
models. And, as in the case of WASP-14b, the remain- 
ing three Spitzer data at 3.6 /fm, 4.5 /tm, and 5.8 are 
also fit equally well by both models, although the 3.6 /tm 
flux is better explained by the C-rich model. The opacity 
sources in the different IRAC channels are the same as 
those discussed for WASP-14b. The photometric data at 
1.6 /tm and 2.1 /tm do not particularly discriminate be- 
tween the two models because both models have similar 
lower atmospheric temperatures, since neither of them 
hosts a thermal inversion. 

However, an important discrimination is provided by 
the z'-band observation (Burton et al. 2012) centered at 
0.9 /xm, in which TiO is a prominent opacity source for 
the 0-rich model. Even though TiO and VO may not be 
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Fig. 13. — Observations and model spectra of the dayside atmo- 
sphere of hot Jupiter WASP-33b. The observed planet-star flux 
contrasts are shown in blue and include photometric data from 
ground-based observations at 0.9 fim (Smith et al. 2011) and at 2.1 
fim (Deming ct al. 2012), and space-based Spitzer observations in 
the IRAC bandpasscs at 3.6 fim and 4.5 ^m (Deming et al. 2012). 
Two model spectra with different C/O ratios and temperature pro- 
files are shown in the solid curves. The data are best explained by 
a carbon-rich model without a thermal inversion, shown in red. 
However, an oxygen-rich model with a thermal inversion can also 
explain all but the 2.1 /im point. The dotted grey curves show 
blackbody spectra of the planet at 2000 K, 2800 K, and 3800 K. 

abundant enough in the upper-atmosphere, due to either 
gravitational settling or chromospheric activity, in an O- 
rich atmosphere it should still be present in the lower 
atmosphere around P ^ 1 bar where the temperatures 
arc hottest and where the incident irradiation doesn't 
reach with significant intensity. Consequently, if C/0 j 

1 and there is no thermal inversion, highly irradiated at- 
mospheres like WASP-19b should display noticeable TiO 
and VO absorption between 0.7 and ^ l.lmicron, which 
is observable as decreased flux in the z', z, and J bands, 
as can be seen in Fig. [T^l The flux observed in the z' 
band is better fit by the C-rich model which doesn't have 
significant TiO /VO absorption. Consequently, we tenta- 
tively classify WASP-19b as a C2 class planet in Fig.[5]in 
view of the marginally better fit provided by the C-rich 
model at 0.9 /tm and 3.6 /tm. 

New observations arc required to conclusively deter- 
mine the C/0 ratio of WASP-19b and to classify it as 

02 or C2. As shown in Fig. [121 new ground-based data 
in the z and J bands will be able to discriminate between 
the different ffux continua at 1 /im and 1.2 /im predicted 
by the 0-rich and C-rich models. Additionally, observa- 
tions with HST WFC3 would be able to provide a robust 
discrimination between the two models in the H2O bands 
between 1.1 /tm and 1.7 /tm. 

5.5. WASP-SSb 

The transiting hot Jupiter WASP-33b, orbiting an A5 
main-sequence star at an orbital separation of 0.03 AU 
is one of the most highly irradiated hot Jupiters known ( 
Collier- Cameron et al. 2010; Herero et al. 2011). At an 
incident irradiation exceeding 10^° ergs/cm^/s, and an 
equilibrium temperature exceeding 3300 K, it the hottest 
planet in our sample and is comparable in temperature 
to low mass stars. Consequently, WASP-33b is the most 
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suitable hot Jupiter to study thermal and chemical char- 
acteristics of irradiated giant planets. If the atmosphere 
of WASP-33b is oxygen-rich, TiO and VO should be 
abundant in the atmosphere, despite any gravitational 
settling. Consequently, an oxygen-rich WASP-33b would 
offer the most unambiguous opportunity to detect a ther- 
mal inversion in a hot Jupiter atmosphere and place it in 
the 02 class discussed in section |4l On the other hand, 
the lack of a thermal inversion in WASP-33b would be 
a strong indication of a carbon-rich atmosphere, i.e a 
C2-class planet. 

Observations of the dayside atmosphere of WASP-33b 
nominally allow both possibilities, 02 as well as C2, thus 
necessitating new observations to break the degeneracy. 
The available data comprise of ground-based photometry 
in the z' band (Smith et al. 2012) and the Ks band (Dom- 
ing et al. 2012), and Spitzer IRAC photometry in the 3.6 
fim and 4.5 fim channels (Deming et al. 2012). Deming 
et al. (2012) conclude that all the existing observations 
are consistent with either a carbon-rich atmosphere with- 
out a thermal inversion (i.e. a C2-type atmosphere) or an 
oxygen-rich atmosphere with a thermal inversion (i.e. an 
02-type atmosphere). While our results in the present 
work generally agree with the assessment of Deming et 
al., we note that a C2-type atmosphere for WASP-33b 
explains the data marginally better than one of 02 type. 
The data and two model spectra of WASP-33b are shown 
in Fig. [131 Table [T] shows the chemical compositions of 
the models. 

The data in the z-band and in the two Spitzer bands, at 
3.6 /im and 4.5 /xm, can be explained almost equally well 
by both models (C2 and 02), but for very different rea- 
sons. The high flux observed in the z' band is naturally 
explained by the higher lower atmospheric temperature 
in the non-inversion model. And since that model is C- 
rich, no significant TiO or VO absorption is expected 
and hence the high predicted flux in the z' band. On 
the other hand, the 02 model which has a thermal in- 
version has a cooler lower lower atmosphere than the 
non-inverted model. However, because this model is O- 
rich substantial TiO and VO is present and is seen as an 
emission feature, as opposed to absorption, due to the 
thermal inversion. Consequently, the 02 model predicts 
a higher brightness temperature in the z' band than its 
lower atmospheric temperature and hence explains the 
high observed z' band flux. Similar reasoning explains 
the fits of both models to the Spitzer data, at 3.6 /im and 
4.5 /xm. At 3.6 /xm the low observed flux is explained by 
the inversion model based on the lower temperature of 
its lower atmosphere and lack of any signiflcant opacity 
source at 3.6 /im in its oxygen-rich atmosphere for T > 
2500 K. The same model explains the higher 4.5 /im flux 
due to CO in emission. On the other hand, the non- 
inverted C-rich model explains the 3.6 /im point due to 
strong absorption by HCN, CH4 and C2H2 in the 3.6 /xm 
channel and absorption due to CO and HCN in the 4.5 
/.jm channel. 

The Kj-band flux, however, favors a C-rich atmosphere 
without a thermal inversion. The 02 model predicts over 
2-a lower planet-star flux contrast than that observed, 
whereas a C2 model fits the data point rather precisely, 
within the 1-a errors. The Kg band is devoid of strong 
molecular absorption and hence probes thermal emission 
from the deepest layer of a planetary atmosphere (P ~ 1 



bar) until continuum opacity makes the atmosphere op- 
tically thick. Consequently, the difference between the 
Kg band fiuxes predicted by the two models directly re- 
flects the difference between the isothermal temperatures 
of the deep atmosphere in the two models. More gener- 
ally, since temperature profiles with thermal inversions 
have cooler lower atmospheres compared to those with- 
out thermal inversions, for the same irradiation, thermal 
emission in opacity windows such as the J, H, and K 
bands tends to be typically lower for models with ther- 
mal inversion. 

WASP-33b is one of the most ideal targets for follow-up 
observations with existing instruments. The very high 
planet-star flux contrast implies that thermal emission 
from the planet can easily be detected with routine obser- 
vations in the ncar-IR as has already been demonstrated 
by Smith et al and Deming ct al. New ground-based ob- 
servations in the J-band (1.2 /im) and H-band (1.6 /tm) 
and 77S'T observations with the WFC3 instrument should 
be able to conclusively distinguish between the two pos- 
sibilities of an 02 and a C2 atmosphere. As shown in 
Fig. [131 J a-iid H band observations should be able to 
place tight constraints on the temperature in the lower 
atmosphere and hence the presence of a thermal inver- 
sion, given the large difference in flux contrast between 
the two models in those bands. Furthermore, an oxygen- 
rich WASP-33b must show strong emission features of 
H2O in the HST WFC3 bandpass due to a thermal in- 
version in accordance with the TiO/VO hypothesis of 
Fortney et al. (2008). On the other hand, weak features 
in the water bands would indicate a carbon-rich atmo- 
sphere in WASP-33b. 

5.6. WASP-12b 

The atmosphere of the transiting hot Jupiter WASP- 
12b (Hebb et al. 2009) has generated substantial interest 
in the recent past. Being one of the most highly irradi- 
ated hot Jupiters, with an equilibrium temperature of 
~ 2600 K, initial models assuming solar abundance (i.e. 
0-rich) compositions predicted strong thermal inversions 
to be present in the dayside atmospheres of planets like 
WASP-12b (Fortney et al. 2008). Even considering de- 
pletion due to gravitational settling, TiO was expected 
to be abundant in the dayside atmosphere of WASP- 
12b and, hence, cause a thermal inversion (Spiegel et al. 
2009). Furthermore, an 0-rich atmosphere would imply 
the presence of abundant H2O and negligible CH4 (Bur- 
rows & Sharp 1999). However, observations of the day- 
side atmosphere of WASP- 12b suggested a stark contrast 
to previous theoretical predictions which assumed solar 
abundances. Using multi-band photometric observations 
over a wide wavelength range (1-10 /im) obtained with 
Spitzer (Campo et al. 2011) and from ground (CroU et 
al. 2010), Madhusudhan et al. (2011a) reported the de- 
tection of a carbon- rich atmosphere (C/0 > 1), implied 
by an under- abundance of H2O and over-abundance of 
CH4, and the lack of a strong thermal inversion in the 
dayside atmosphere of WASP-12b. Consequently, it is 
apparent that WASP-12b cannot be characterized solely 
on the basis of incident irradiation, and assuming a so- 
lar C/0 ratio, as per previous classification schemes (e.g. 
Fortney et al. 2008). However, the C/0 > 1 and the lack 
of a strong thermal inversion in WASP- 12b are both con- 
sistent with the properties of a C2-class planet according 
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Fig. 14. — Observations and model spectra of dayside thermal 
emission from the hot Jupiter WASP-12b. The data are shown 
in blue and black circles with error bars. The blue data show 
ground-based broadband photometric points at 0.9 fira (Lopez- 
Morales et al. 2010), 1.2 ^m, 1.6 fim, and 2.1 fim (CroU et al. 
2010), and Spitzer IRAC photometry at 3.6 fim, 4.5 fim, 5.8 ^m 
and 8 ^m (Campo ct al. 2011, Cowan et al. 2012). The data 
in black circles show a near-infrared spectrum obtained using the 
HST Wide Field Camera 3 (Swain et al. 2012). The photometric 
bandpasscs are shown in black solid curves at the bottom of the 
plot and the HST WFC3 bandpass is shown in the inset. The 
solid curves show two models with different C/O ratios. The data 
are explained well by a carbon-rich model, whereas a model with 
oxygen-rich composition is ruled out by the data, consistent with 
the findings of Madhusudhan et al. (2011a). Both models have 
the same temperature profile without a thermal inversion; a strong 
thermal inversion is ruled out by the data. The dotted gray curves 
show the planet blackbody spectra at 1800 K. 2500 K. and 3200 
K. 

to our 2-D classification scheme suggested in section 

New data confirm a high C/O ratio and a weak ther- 
mal inversion in the dayside atmosphere of WASP- 12b. 
Using new observations in the Spitzer IRAC channels at 
3.6 ^m and 4.5 fxra, in addition to existing data, Cowan 
et al. (2012) confirmed that their best fitting models re- 
quired a high C/O ratio and the lack of a strong thermal 
inversion in the dayside atmosphere of WASP-12b (see 
Section 5.2.1 of Cowan et al.). It is to be noted, how- 
ever, that since Cowan et al. used forward models (Bur- 
rows et al. 2008) for their interpretation, as opposed to 
a retrieval algorithm of Madhusudhan et al, a statistical 
estimate of the C/O ratio was not possible in the Cowan 
et al study, despite their confirmation of a C/O > 1 im- 
plied by their high CO/II2O requirement. New studies 
have also reported ground-based photometric (Zhao et al. 
2012) and spectroscopic (Crossfield) observations of ther- 
mal emission from WASP-12b in the K band (around 2.1 
^m), and confirmed previous measurements in the same 
bandpass by CroU et al. (2010). More recently. Swain 
et al. (2012) reported observations of WASP-12b in the 
near- infrared (1-1.7 fj,m) using the HST Wide Field Cam- 
era 3 (WFC3) . Based on their thermal emission spectra, 
Swain et al. suggested that even though models with 
C/O > 1 best explained their data, they could not pro- 
vide conclusive constraints on the C/O ratio. However, 
Swain et al. stress that explaining their data required a 
substantial paucity of H2O (< 10"®). It is well known 
that such low levels of H2O is possible if either C/O > 
1, or if the overall metallicity of the planet is low (e.g. 



see Section 3; also see Seager et al. 2005, Madhusud- 
han et al. 2011, Moses et al. 2012). In the latter case, 
both C/H and O/H would have to be under-abundant 
in WASP- 12b by about two orders of magnitude below 
the stellar values, a scenario which would not be able 
to explain the high CO required to explain the Spitzer 
data in Madhusudhan et al. (2011a) and Cowan et al. 
(2012), which Swain et al did not consider including in 
their model fits. Furthermore, a hundred times lower 
metallicity (particularly C and O abundances) in a giant 
planetary atmosphere relative to its host star would be 
challenging to explain using the standard core accretion 
model for giant planet formation (see e.g. Owen et al. 
1999; Atreya & Wong 2005). 

In the present work, we find that all existing observa- 
tions of thermal emission from WASP- 12b are consistent 
with a high C/O ratio and a weak thermal inversion in 
its dayside atmosphere, confirming the findings of Mad- 
husudhan et al. (2011a). The data include four chan- 
nels of ground-based near-infrared photometry in the z' 
(Lopez-Morales et al. 2010), J, H, and K bands (Croll 
et al. 2011), four-channels of Spitzer IRAC photome- 
try (Campo et al. 2011; for 3.6 /im and 4.5 /im wc use 
the new values from Cowan et al. 2012), and a near- 
infrared spectrum in the HST WFC3 bandpass (Swain 
et al. 2012). We have also updated our spectral models, 
over those used in Madhusudhan et al. (2011), by includ- 
ing molecular opacities due to C2II2 and HCN, which 
could be important sources of opacity in carbon-rich at- 
mospheres (see Section 3). Fig. [HI shows all the data 
along with two spectral models, and the chemical com- 
positions of the models are shown in Table [TJ We find 
that a carbon-rich model (with C/O > 1) provides the 
best fit to all the data whereas a solar abundance model 
is inconsistent with the data, in agreement with the in- 
terpretation of Madhusudhan et al. (2011a). The pri- 
mary difference between the fits in Madhusudhan et al. 
(2011a) and the present work lies in the interpretation of 
the absorption in the 8^m IRAC channel. While Mad- 
husudhan et al. (2011a) suggested that the low flux at 
8/im was solely due to CII4 absorption, we find that the 
low flux could be due to a combination of CII4 and HCN 
opacity. The uncertainties in high temperature opacities 
of hydrocarbons (Tinetti et al. 2010) makes it difficult to 
judge the exact relative contributions of CH4 and HCN, 
but it is possible that in C-rich atmospheres HCN may 
dominate the 8^m opacity, as has also been suggested by 
Kopparapu et al. (2012) and Moses et al. (2012). 

We classify WASP-12b as a C2 class planet, based on 
all observations of its dayside atmosphere. The HST 
WFC3 spectrum reported by Swain et al. (2012) pro- 
vides a particularly unambiguous constraint on the C/O 
ratio as shown in the inset of Fig. [TH The WFC3 data 
by themselves alone conclusively rule out H2O absorp- 
tion (as also reported by Swain et al. 2012) which is a 
prominent feature around 1.4 ^m, and hence rule out an 
0-rich atmosphere. Furthermore, the small absorption 
feature in the data at 1.57 ^m appears to be consis- 
tent with the small C2H2 feature in the C-rich model, 
although the detection is only marginal. The data also 
rule out a strong thermal inversion in the atmosphere, 
which is apparent from the lack of any strong emission 
features in the spectrum, as shown in Fig. [M] 

Additional constraints on the atmospheric composition 
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of WASP- 12b are expected in future work. In this work, 
we have not considered observations of transmission spec- 
tra of WASP- 12b which probe the atmospheric compo- 
sition at the day-night terminator of the planet. Cowan 
et al. (2012) reported transmission photometry using 
Spitzer at 3.6 /im and 4.5 /im. They observed unusu- 
ahy high absorption in the 3.6 /xm channel compared to 
the the 4.5 /.tm channel, which they were unable to ex- 
plain using any composition in their atmospheric mod- 
els, in which the dominant sources of opacity were only 
H2O and CO, and, probably, CH4. Swain et al. (2012) 
reported transmission spectroscopy of WASP- 12b using 
HST WFC3. In order to explain their data. Swain et al. 
required extremely low levels of H2O (^ 5 x 10~^) and 
a possible presence of TiH and CrH in the atmosphere. 
Such a composition is possible if either the atmospheric 
metallicity is at least two orders of magnitude lower than 
the stellar metallicity or if the metallicity is normal but 
the C/0 ratio is > 1 as has been inferred on the day- 
side by Madhusudhan et al. (2011) and in the present 
work. Future observations of WASP- 12b in transit in 
the near- infrared (between 0.7 fim and 1.1 /im) could 
also place constraints on the TiO/VO abundances in its 
atmosphere. In an upcoming effort, we plan to reanalyze 
all existing transmission data of WASP-12b and fit them 
with our models which include a broader set of opacities, 
than previous studies, including HCN and C2H2; to ob- 
tain statistical constraints on the atmospheric C/0 ratio 
at the day-night terminator of WASP-12b. 

6. OBSERVING STRATEGY AND OPTIMAL SYSTEMS 

Future observations in the near-infrared with existing 
facilities will be able to constrain C/0 ratios for sev- 
eral currently known exoplanets. As discussed for each 
system in Section [SJ robust constraints on the C/0 ra- 
tios and thermal structures of hot Jupiter atmospheres 
can be placed currently from thermal spectra using a 
combination of existing Spitzer and ground-based data 
with new data from HST WFC3 and additional ground- 
based observations. The HST WFC3 bandpass would be 
critical to constrain the H2O abundances in hot Jupiter 
atmospheres, especially for planets for which a Spitzer 
5.8 fim observation is not available. At the same time, 
ground-based near-infrared observations provide a very 
powerful means to constrain the temperature structure 
independent of the composition. The J, H, and K bands 
probe deeper in the planetary atmosphere, compared to 
the mid-IR Spitzer bands, and hence constrain the tem- 
perature structure in the lower atmosphere, which is im- 
portant to constrain a thermal inversion in the upper at- 
mosphere. Atmospheres with thermal inversions in the 
upper-atmospheres tend to have lower temperatures in 
the lower atmospheres, and vice versa (see e.g. Bur- 
rows et al. 2008, Madhusudhan & Seager, 2010). Once 
the temperature profile is partially constrained by the 
near-infrared data, the Spitzer data along with any ad- 
ditional observations in the molecular bands, allow good 
constraints on the molecular composition and hence the 
C/0 ratio. 

For very hot atmospheres, the near-infrared region be- 
tween 0.8 and 1.1 fxui also provide an important diag- 
nostic for the temperature profile as well as the C/0 
ratio. This wavelength region contains strong spectral 
features due to TiO and VO. As discussed in the case of 
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Fig. 15. — Transiting hot Jupitcrs that arc optimal for charac- 
terizing atmospheric C/O ratios (see Section [6J. The left vertical 
axis shows the incident flux from the host star received by the 
planet at the sub-stellar point. The right vertical axis shows the 
equilibrium temperature with zero albedo and no redistribution 
(Toq f— 0), which is a nominal lower limit on the quench tempera- 
ture (Tq), as discussed in Section 13.31 All the 21 systems shown 
here have temperatures greater than 1300 K in their lower atmo- 
spheres, at P ~ 1 bar, and have all been observed in at least three 
channels of Spitzer photometry. Combining existing data of these 
systems with new observations in the near-infrared, that are pos- 
sible with existing facilities, can yield constraints on their C/O 
ratios. The systems parameters for the planets are adopted from 
|http: / /www exoplanets. org| 

WASP-19b, in Section O the z' band, Y band, z band, 
and part of the J band, all between 0.9 and 1.1 ^m, can 
be useful probes of the TiO abundance and, indirectly 
the C/0 ratio. If a highly irradiated atmosphere, such 
as WASP-19b, is known to have no thermal inversion, 
based on longer wavelength data, then one or more of 
these four bands {z' to J) can be used to distinguish 
between a C-rich and an 0-rich atmosphere. A highly ir- 
radiated 0-rich atmosphere without a thermal inversion 
will have strong TiO absorption, and hence low flux, in 
these bands. Whereas a C-rich planet without a thermal 
inversions can have high flux in these bands, depending 
on how much TiO and VO are depleted due to the high 
C/0. Observations in the same bands for a model with a 
thermal inversion are discussed for the case of WASP-33b 
in Section [SJ 

New observations aiming to determine C/0 ratios will 
benefit most from observing the hottest of hot Jupitcrs. 
As discussed in Section 3, the contrast in molecular com- 
positions between C-rich and 0-rich atmospheres occurs 
increases with temperatures. While the contrast is min- 
imal for temperature below ~ 1200i4r, the contrast is 
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maximal for temperatures above ~ 2000 K. Direct ele- 
mental estimates are possible if observations can be made 
in the H2O, CO, and/or CH4 bands, and, if possible, in 
bands of C2H2 and HCN. The inference of C/0 ratios 
is straightforward at high temperatures. A strong H2O 
feature and/or a weak CH4 feature is an indication of 
C/0 < 1. Conversely, a weak H2O feature and/or a 
noticeable CH4 feature is an indication of C/0 > 1. Ob- 
servations in CO bands should always find a strong CO 
signature for hot Jupiters, independent of the C/0 ra- 
tio. The CO/H2O ratio wiU be unity for C/O=0.5, and 
can be as high as 10-1000 for higher C/0 ratios. More 
detailed inferences based on model spectra arc discussed 
in sections 13.41 and [5l 

We identify a sample of 20 hot Jupiters which are cur- 
rently most suitable for measuring C/0 ratios, as shown 
in FigHSl These planets all have Tqj^o si 1300, and have 
already been observed in at least three channels of Spitzer 
photometry, in the cryogenic phas^H The Spitzer channels 
cover features of all the prominent molecules: 3.6 //m 
(H2O and CH4, C2H2, HCN), 4.5 Aim (CO, CO2, C2H2, 
HCN), 5.8 Aim (CO and H2O), 8 /im (H2O, CH4, C2H2, 
HCN), 16 /im (CO2 and H2O), and 24 //m (H2O). As 
such, combining existing Spitzer data for systems shown 
in FigllSI with new ncar-IR observations with HST and 
ground-based facilities can lead to C/0 estimates for a 
statistically significant sample of 20 hot Jupiter atmo- 
spheres. 

High confidence estimates of C/0 ratios can be made 
either with spcctrophotomctric observations that can re- 
solve molecular bands or by combining several broad- 
band photometric observations over a wide wavelength 
range, or both. Currently available spectrophotomet- 
ric facilities include the HST WFC3 instrument in the 
1 — 1.7 /im bandpass (Doming et al. 2010 HST Program 
#12736), ground-based spectrophotometry (e.g. Swain 
et al. 2010; Mandell et al. 2011; Crossficld et al. 2011; 
Beanetal. 2011), and SOFIA (Angerhauscn ct al. 2010), 
in the near term, and the James Webb Space Telescope 
in the future. On the other hand, photometric obser- 
vations can be made with Warm Spitzer at 3.6 /im and 
4.5 /xm, and ground-based ncar-IR photometry in the z' 
(Lopez-Morales et al. 2010), J, H, and Ks (CroU et al. 
2011). 

While the above discussion focused primarily on ob- 
servations of thermal emission from exoplanets, trans- 
mission spectra can also be used to constrain the C/0 
ratios in giant exoplanetary atmospheres. Since trans- 
mission spectra arc less sensitive to temperature struc- 
ture, and hence to thermal inversions, the benefits of 
observing in continuum such as broadband photometry 
in H or K bands is minimal, except to measure a refer- 
ence absorption or an accurate photosphcric radius of the 
planet. However, transmission spectroscopy or narrow- 
band photometry in the molecular bands would be a 
powerful means to constrain the chemical composition. 
Some of the useful bands in that regard would be the 
TJST WFC3 bandpass between 1.1 and 1.7 /im for H2O, 
the 0.8 - 1.1 Aim region for TiO/VO, parts of the H (1.6 
/im) and K (2.1 /im) bands for CH4 and other hydro- 
carbon features, and the CO band around ^ 2.3 /.tm. If 

^ In the post-cryogenic "warm" phase now only two of Spitzer 
's original six photometric channels are operational. 



either H2O or CH4 abundance can be measured in trans- 
mission, the C/0 ratio can be estimated. If only CO is 
detected, it would not be adequate to determine the CO 
ratio. Measurement of any hydrocarbon abundance in a 
hot atmosphere would be a strong evidence of a carbon- 
rich atmosphere. 

7. SUMMARY AND CONCLUSIONS 

We propose the C/0 ratio as a new dimension for at- 
mospheric characterization of extrasolar planets. In the 
past, atmospheric models for exoplanets have typically 
assumed solar elemental abundances. Such models have 
been used for one-dimensional characterization of irradi- 
ated giant planets based primarily on the degree of inci- 
dent irradiation, or equivalently temperature, while fix- 
ing the elemental abundances to solar values, with C/0 
= 0.5. However, recent observations of cxoplanctary at- 
mospheres have revealed extreme departures from pre- 
dictions based on models and classification schemes that 
assumed solar abundances. In this work, we find that 
the C/0 ratio plays a critical role in governing the atmo- 
spheric chemistry, and non-solar C/0 ratios can poten- 
tially explain observations which were hitherto reported 
as anomalous using solar abundance models. We, there- 
fore, propose a two-dimensional classification system for 
cxoplanctary atmospheres where chemistry, represented 
by the C/0 ratio, and temperature (T), or equivalently 
incident irradiation, form the two dimensions. 

Our proposed 2-D characterization scheme opens a new 
physico-chemical phase space in which to interpret obser- 
vations of H2-dominated exoplanetary atmospheres. In 
hot H2-dominated atmospheres, the C/0 ratio critically 
infiuences the concentrations of major molecular species, 
such as H2O, CO, and CH4. For T > 1200 K, a C/0 
of 1 constitutes a natural boundary between the famil- 
iar 0-rich (C/0 < 1) and the exotic C-rich (C/0 > 1) 
atmospheres. Hot 0-rich atmospheres are expected to 
be dominated by H2O and CO, and be less abundant in 
CH4; CH4 decreases as T increases and becomes negli- 
gible for T > 1500 K. On the other hand, C-rich atmo- 
spheres at high T are expected to be depleted in H2O, 
and enhanced in CH4, C2H2, HCN, and in the CO/H2O 
ratio, compared to those obtained with a solar C/0 of 
0.5. Each of the 0-rich and C-rich regimes can be further 
divided into two subclasses depending on the tempera- 
ture or incident irradiation, analogous to previous classi- 
fication schemes (Fortney et al. 2008). However, C-rich 
atmospheres are unlikely to host thermal inversions due 
to TiO/VO even for the most highly irradiated cases, as 
demonstrated by Madhusudhan et al. (2011b). For C/0 
> 1, almost all the oxygen is occupied by CO, which is 
the dominant O bearing molecule at high T, leaving little 
O for other oxygen-bearing molecules. 

We identify four classes of irradiated atmospheres (01, 
02, CI, and C2) with distinct compositional and spec- 
troscopic properties. Classes 01 and 02 are 0-rich (C/0 
< 1), and classes CI and C2 are C-rich (C/0 > 1). The 
division between 01 and 02, and between CI and C2 is 
based on irradiation (or temperature); 01 and CI are the 
cooler classes. In 01 and 02 atmospheres, H2O is a ma- 
jor oxygen-carrier at all temperatures, and CO or CH4 
is the dominant carbon-carrier, depending on the tem- 
perature. The division between the 01 and 02 classes 
is analogous to the pL and pM classes, respectively, of 
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Fortnoy et al. (2008). 02 classes might be likely to host 
thermal inversions due to TiO because of their high ir- 
radiation levels. On the other hand, CI and C2 classes 
are characterized by depleted H2O and enhanced CH4. 
C2H2; and HCN compared to 0-rich atmospheres at the 
same level of irradiation. The differences in H2O and 
CH4 abundances between 0-rich and C-rich atmospheres 
are enhanced with temperature and can reach upto fac- 
tors of 10^ for T > 2000 K. In the cooler CI class, 
however, the H2O depletion is too low to be distinguish- 
able from atmospheres with the same irradiation in the 
01 class. Consequently, C/0 determination of CI class 
atmospheres would require high resolution spectroscopy 
and/or other molecular species as tracers of C/0 > 1. 
Finally, neither CI nor C2 class atmospheres can host 
thermal inversions due to TiO or VO. 

We emphasize the importance of considering non-solar 
C/0 ratios in interpreting observations of exoplanetary 
atmospheres. We find that non-solar C/0 ratios can 
potentially explain observations of several hot Jupiters 
which have hitherto been regarded as anomalous in the 
literature, based on solar abundance models. For exam- 
ple, while observations of one of the coolest hot Jupiters, 
XO-lb, were interpreted as indication of a thermal in- 
version (Machalek et al. 2008), some of the hottest 
hot Jupiters like WASP-12b, WASP-14b, and WASP- 
19b showed no evidence of a strong thermal inversion 
(Madhusudhan ct al. 2011a). On the other hand, pre- 
vious studies have been unable to explain observations 
of CoRoT-2b. We find that the data for both XO-lb 
and CoRoT-2b, along with several other hot Jupiters 
can be explained by models with non-solar C/0 rations 
and with different constraints on their temperature struc- 
tures. These solutions agree with predictions based on 
the 2-D characterization scheme proposed in this work. 
The candidate C-rich atmospheres of XO-lb, CoRoT-2b, 
WASP-19b, and WASP-33b, reported in this work can 
be confirmed with follow-up observations with existing 
facilities. 

We present observational strategics to constrain atmo- 
spheric C/0 ratios of currently known exoplancts with 
existing and forthcoming facilities. The differences be- 
tween C-rich and 0-rich systems increase with temper- 
ature. As such, the hotter atmospheres are more fa- 
vorable for obtaining C/0 constraints. We identify 20 
currently known transiting giant exoplancts which are 
most suitable for C/0 constraints and 2-D character- 
ization. The identified planets all have T > 1200 K in 
their lower atmospheres, and have already been observed 
in at least three channels of Spitzer IRAC photometry. 
As demonstrated for the case of WASP- 12b (Madhusud- 
han et al. 2011a), combining the Spitzer observations 
with near-infrared observations can provide statistically 



meaningful joint constraints on C/0 ratios and temper- 
ature profiles. Near-infrared observations may include 
ground-based photometry (e.g. in the z, J, H, and K 
bands; Lopez-Morales et al. 2010; Croll et al. 2011) and 
spectroscopy (Bean et al. 2010; Mandell et al. 2010; 
Swain et al. 2010; Crossfiled et al. 2011), and space- 
borne spectroscopy with HST WFC3 (e.g. Deming et 
al. 2010 HST Program #12736) and, potentially, SOFIA 
(Angerhausen et al. 2010). More stringent constraints 
will be possible with the James Webb Space Telescope in 
the longer term. 

If atmospheric elemental abundances reflect the abun- 
dances in planetary interiors, exoplancts in the 0-rich 
and C-rich regimes are expected to span a diverse range 
of interior compositions and planet formation scenar- 
ios. In particular, C-rich atmospheres initiate the new 
class of carbon-rich planets (CRPs; Madhusudhan ct al. 
2011b) with atmospheres, interiors, and formation sce- 
narios markedly different from expectations based on the 
solar abundances. While 0-rich chemistry is expected 
to form silicates in planetesimals and rocky cores, C- 
rich environments would be expected to be dominated 
by carbides (such as SiC) and pure carbon minerals such 
as graphite and diamond (Gilman 1969; Larimer 1974; 
Kuchncr & Seager, 2005; Bond et al. 2010). Equally 
drastic are the protoplanetary conditions required to 
form CRPs, such as a naturally carbon-rich gas in the 
disk (Obcrg et al. 2011), or a sink of oxygen in the disk 
(Madhusudhan et al. 2011b), or the predominance of 
planetesimals rich in carbonaceous material like tar, in- 
stead of water-ice (Lodders, 2004). Thus, determination 
of C/0 ratios in exoplanetary atmospheres are critical 
for understanding their interiors and formation scenar- 
ios, which have thus far been commonly studied assum- 
ing solar abundances. 
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